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Abstract 


Free  electron  laser  theory  is  developed  from  the  Maxwell  and  Lorentz  force  equations; 
the  properties  and  characteristics  of  the  laser  are  reviewed.  The  wave  equation  is  solved 
for  the  fundamental  Gaussian  mode,  and  higher-order  modes  in  Cartesian  and  cylindrical 
coordinate  spaces,  yielding  expressions  for  the  complete  and  orthogonal  basis  sets  of  Hermite- 
and  Laguerre-Gaussian  beams.  Motivated  by  the  evident  inclusion  of  higher-order  modes 
in  free  electron  laser  simulations,  a  tool  is  developed  for  the  higher-order  (in  particular 
Laguerre-Gaussian)  modal  analysis  of  simulated  free  electron  laser  beams. 


v 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


vi 


Contents 


I.  The  Free  Electron  Laser  1 

A.  System .  1 

B.  General  Theory .  3 

1.  Electron  Motion .  3 

2.  Evolution  of  the  Optical  Mode .  9 

3.  Interaction  and  Gain .  12 

II.  Fundamental  Mode  Optics  19 

A.  Propagation  -  Analytic  Description .  19 

1.  The  Parabolic  Wave  Equation  .  19 

2.  Solving  the  Dimensionless  Parabolic  Wave  Equation  ....  21 

B.  Properties  of  the  Gaussian  Beam  .  24 

IIIMOTIVATION:  FREE  ELECTRON  LASER  SIMULATIONS  27 

A.  Multiple  Modes .  27 

1.  Multi-Mode  Simulations  of  a  Short-Rayleigh  Length  FEL  27 

2.  Results .  28 

B.  Optical  Mode  Distortion .  34 

1.  Optical  Mode  Distortion  in  a  Short  Rayleigh  Length  FEL  34 

2.  Results .  36 

C.  The  Need  for  Higher- Order  Mode  Analysis .  37 

IV .Higher-Order  Modes  39 

A.  Hermite-Gaussian  Beams  .  39 

B.  Laguerre- Gaussian  Beams .  41 

C.  Numerical  Propagation  of  Higher-Order  Modes .  47 

D.  Power  on  Target .  49 

E.  Conclusion .  53 

Appendix:  The  Parabolic  Wave  Equation  55 


Initial  Distribution  List 


59 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


List  of  Figures 


1.1  Oscillator  FEL  .  1 

1.2  Amplifier  FEL .  2 

1.3  Recirculating  free  electron  laser  oscillator .  2 

1.4  FEL  Gain  on  Resonance .  8 

1.5  FEL  Gain  Off  Resonance .  8 

1.6  FEL  Weak  Field  Gain  Curve .  16 

2.1  Optical  Mode  Diagram .  26 

2.2  Propagation  of  a  Gaussian  beam .  26 

3.1  Output  from  a  2003  FEL  simulation .  29 

3.2  Extraction  77  vs.  bunch  charge  q .  30 

3.3  Extraction  rj  vs.  undulator  periods  N .  31 

3.4  Extraction  77  and  optimum  initial  phase  velocity  uq  vs.  normalized  electron 

beam  waist  radius  a .  32 

3.5  The  interaction  region .  32 

3.6  Extraction  77  and  mirror  intensity  vs.  normalized  Rayleigh  length  50 .  33 

3.7  Output  from  a  2004  FEL  simulation .  35 

3.8  Extraction  77  vs.  dimensionless  Rayleigh  length  Zo .  36 

3.9  Extraction  77  vs.  dimensionless  mirror  separation  rm -Ir .  36 

3.10  Optical  mode  shape  suggesting  reinforcement  by  electron  beam .  37 

4.1  Hermite-Gaussian  Modes .  41 

4.2  Laguerre- Gaussian  Modes  .  43 

4.3  Output  from  bm.c .  48 

4.4  Output  from  bm.c .  49 

4.5  Two  functions  lead  to  the  same  squared  function .  49 

4.6  Output  from  bm.c .  50 

4.7  A  combination  of  LG  modes  before  and  after  propagation .  50 

4.8  Output  from  bm.c .  51 

4.9  Output  from  bm.c .  51 

4.10  Energy  on  target  for  higher-order  modes .  51 

4.11  Three-dimensional  plot  of  L\ .  52 

4.12  Various  beam  modes  propagating .  52 


IX 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


x 


List  of  Tables 


2.1  Symbols  Used  in  the  Dimensionless  Parabolic  Wave  Equation .  25 

3.1  Parameters  for  a  Proposed  High- Power  FEL .  34 

4.1  Hermite  Polynomials .  41 

4.2  Laguerre  Polynomials .  42 


THIS  PAGE  INTENTIONALLY  LEFT  BLANK 


List  of  Symbols 


a 

A 

Ab 

A0 

B 

Br 

Bu 

c 

e 

E 

Eb 

Ee 

Er 

F 

paj3 

G 

J 

J 

k 

k0 

K 

h 

L 

m 

N 

Ne 

q 

Q 

rb 

t 

u 

Ua 

V 


dimensionless  optical  field  strength 
vector  potential 

electron  beam  cross  secitonal  area 
optical  mode  cross  sectional  area 
four-vector  potential 
magnetic  field 

optical  mode  magnetic  field  magnitude 

undulator  magnetic  field  magnitude 

speed  of  light  in  vacuum 

electron  charge 

electric  held 

electron  beam  energy 

electron  energy 

optical  mode  electric  held  magnitude 
filling  factor 
hcld-strength  tensor 
gain 

current  density 
dimensionless  current  density 
four-current  density 
optical  wavenumber 
undulator  wavenumber 
undulator  parameter 
bunch  length 
undulator  length 
electron  mass 

number  of  undulator  periods 
number  of  electrons 
charge 

resonator  quality  factor 

beam  radius 

time 

electron  velocity 
electron  four- velocity 
volume 


xiii 


w  spot  size 

w o  spot  size  at  waist 

zo  Rayleigh  range 

Zq  dimensionless  Rayleigh  range 

1 3  dimensionless  velocity 

7  Lorentz  factor 

£k  polarization  unit  vector 

(  ponderomotive  electron  phase 

rj  extraction 

A  optical  wavelength 

A0  undulator  wavelength 

v  dimensionless  phase  velocity 

a  dimensionless  electron  beam  radius 

pe  electron  density 

r  dimensionless  time 

T  proper  time 

0, 0  optical  phase 

uj  optical  frequency 

cuq  undulator  frequency 


xiv 


Chapter  I. 


The  Free  Electron  Laser 


In  1971  Stanford  University  physicist  John  Madey  proposed  the  development  of  a  coherent 
light  source  based  on  the  radiation  emitted  by  relativistic  free  electrons:  a  free  electron 
laser.1  In  1976,  his  team  amplified  a  CO2  laser  beam  using  relativistic  free  electrons  from 
the  Stanford  Superconducting  Linear  Accelerator.  They  later  observed  coherent  infrared 
radiation  at  a  wavelength  of  3  /irn  in  an  oscillator  configuration.  These  achievements  gave 
birth  to  the  contemporary  free  electron  laser  described  below. 

A.  System 

The  two  main  components  of  a  free  electron  laser  are  an  accelerator  which  imparts  kinetic 
energy  to  a  beam  of  free  electrons,  and  an  undulator  also  called  a  wiggler  which  extracts 
energy  from  the  beam  in  the  form  of  coherent  radiation.  There  are  a  variety  of  FEL  designs; 
all  contain  these  two  essential  elements.  Varying  among  individual  configurations  are  the 
methods  by  which  the  electron  beam  is  created,  transported,  then  dumped,  and  how  the 
radiation  is  amplified.  Some  FEL’s  are  amplifiers,  some  are  oscillators.  The  oscillator 


configuration  is  shown  in  Figure  1.1.  Electromagnetic  radiation  is  stored  within  a  resonant 
cavity,  and  the  increased  optical  power  is  outcoupled  through  a  partially-transmitting  dielec¬ 
tric.  In  the  amplifier  configuration  (Figure  1.2),  spontaneously  emitted  light  or  light  from  a 
seed  laser  is  amplified  in  a  single  pass  through  the  undulating  electron  beam. 

1J.  M.  J.  Madey,  “Stimulated  Emission  of  Brehmsstrahlung  in  a  Periodic  Magnetic  Field”,  Journal  of 
Applied  Physics ,  Vol.  42,  No.  5,  pp.  1906-1913,  April  1971. 
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Electron  Beam 


Figure  1.2:  Amplifier  FEL 


The  source  of  free  electrons  is  usually  a  cathode  within  an  electron  gun.  Whether  con¬ 
figured  as  an  oscillator  or  amplifier,  the  electron  beam  itself  may  be  either  recirculated  or 
dumped  after  exiting  the  undulator.  Recirculation  allows  the  recovery  of  beam  energy  back 
to  the  accelerator,  enhancing  overall  system  efficiency,  with  the  added  benefit  of  dumping 
less  energetic  electrons,  producing  less  background  radiation.  In  the  RF  recovery  configu¬ 
ration,  electrons  return  to  the  accelerator  180°  out  of  phase  with  respect  to  those  entering 
from  the  beam  source.  Figure  1.3  shows  the  entire  system  for  a  recirculating  free  electron 
laser  oscillator. 
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Figure  1.3:  Recirculating  free  electron  laser  oscillator. 


Attributes 

The  free  electron  laser’s  most  striking  difference  from  conventional  lasers  is  also  its  chief 
advantage.  The  absence  of  a  vulnerable  medium  such  as  a  crystal  or  fluid  chemical  matrix 
means  that  an  FEL  can  achieve  extremely  high  intensity  without  damaging  the  laser  itself. 
High  wall-plug  efficiency  (~10%)  is  predicted  for  high  power  FELs  with  recirculating  electron 
beams.  High  reliability  is  also  a  characteristic:  currently  operating  FELs  run  continuously 
and  reliably  for  weeks. 
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Free  electron  lasers  are  continuously  tunable.  The  relativistic  Lorentz  factor  7  of  the 
electron  beam,  the  undulator  period  Ao,  and  undulator  RMS  magnetic  field  strength  £>rms 
determine  the  laser  wavelength  A  according  to  the  resonance  condition2 


A 


A0(l  +  K2) 

272 


where  the  undulator  parameter  K  is  defined  as 


K 


cFArns  A(j 

27rmc2 


(1.1) 


(1.2) 


Here  e  is  the  electron  charge,  m  the  mass  of  the  electron,  and  c  the  speed  of  light  in  vacuum. 
The  relation  (1.1)  suggests  that  the  undulator  wavelength  and  magnetic  field  as  well  as  the 
energy  of  the  electrons3 4 5  may  be  adjusted  within  an  operating  laser  to  yield  light  over  a 
range  of  wavelengths.  Free  electron  lasers  have  been  designed  to  produce  radiation  from 
microwaves  to  X-rays,  and  they  are  continuously  tunable  over  a  smaller  range  within  these 
regimes.  In  addition,  the  short-pulse  nature  of  the  radiation  may  enhance  propagation  of 
optical  power  through  the  atmosphere,  making  the  laser  particularly  suited  as  a  weapon.4,5 


B.  General  Theory 

Free  electron  laser  theory  addresses  the  motion  of  the  oscillating  electrons,  the  propagation 
and  diffraction  of  the  resulting  radiation  (the  optical  mode),  and  the  subsequent  interaction 
of  that  radiation  with  the  electron  beam  itself.  We  forego  a  discussion  of  diffraction  until 
the  next  chapter,  and  proceed  with  a  description  of  the  electron  motion,  the  optical  mode, 
and  gain. 

1.  Electron  Motion 

The  trajectory  of  a  relativistic  charged  particle  in  external  fields  is  governed  by  the  Lorentz 
force  equation6 

dL  =  (1.3) 

clT  me 

2 At  resonance  in  an  FEL,  one  wavelength  of  light  passes  over  one  electron  in  one  undulator  period. 

3Tuning  via  electron  beam  energy  is  generally  considered  more  difficult  than  via  undulator  geometry. 
Indeed,  significant  changes  in  electron  energy  require  adjustments  to  guiding  and  focusing  magnets  along 
the  entire  length  of  the  accelerator. 

4P.  Sprangle,  J.  R.  Penano  and  B.  Hafizi,  “Propagation  of  Intense  Short  Laser  Pulses  in  the  Atmosphere”, 
Phys.  Rev.  E,  66,  046418,  2002. 

5P.  Sprangle,  J.  R.  Penano,  A.  Ting,  B.  Hafizi,  and  D.  F.  Gordon,  “Propagation  of  Short,  High-Intensity 
Laser  Pulses  in  Air”,  J.  Dir.  Energy  1 ,  73,  2003. 

6We  shall  use  T  for  proper  time;  r  will  be  used  for  dimensionless  time;  cgs  units  are  used  throughout  this 
chapter. 
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where  Ua  =  (7 c,  711)  is  the  four- velocity  of  the  particle  and 


paP 


/ 

0 

~EX 

~Ey 

~EZ 

\ 

Ex 

0 

~BZ 

By 

Ey 

Bz 

0 

-Bx 

V 

Ez 

—  By 

Bx 

0 

/ 

(1.4) 


is  the  field-strength  tensor,  whose  components  are  given  by  the  periodic  magnetic  held  of 
the  undulator  and  the  electric  and  magnetic  components  of  the  radiation  held.  For  the 
propagation  of  a  plane  wave  along  the  z-axis  of  a  helical  undulator,7  the  undulator  and 
radiation  helds  are  taken  to  be 


B„ 

=  B0(coskoz,  sinfco^,  0), 

(1.5a) 

Er 

=  FI  (cos  0,  —  sin  0,0), 

(1.5b) 

B, 

=  F?(sin0,  cos0, 0), 

(1.5c) 

where  B0  and  ko  =  27t/Ao  are  the  amplitude  and  wavenumber  of  the  undulator  held,  E  = 
E(x,y,z,t)  is  the  optical  held  amplitude,  and  the  phase  is  0  =  kz  —  ut  +  (j)(z,t ),  where 
k  =  uj/c  is  the  optical  wavenumber,  and  <j>(z,t)  is  the  optical  phase.  The  field-strength 
tensor  is  then 


FaP  = 


(  0  —Ecos'ijj  E  sin  0 

E  cos  0  0  0 

—E  sin  -0  0  0 

V  0  —  Bq  sin  koz—E  cos  ip  Bq  cos  koz+Esinip 


0  \ 

Bq  sin  koz-\-E  cos  ?/> 
—Bo  cos  koz— E  sin  xfj 

0 


Inserting  (1.6)  into  (1.3)  gives  the  four  coupled  differential  equations 
dU° 


dT 
dU 1 
~dT 
dU 2 

hr 

dU 3 
~dT 


= - —  (U1E  cos  ip  —  U2E  sin  ip)  , 

me  x 

= - [ U°E  cos  ip  —  U3  (Bo  sin  koz  +  E  cos  0)1  , 

me  1  J 

=  — —  \—U°E  sin0  —  U 3  (— Bocoskoz  —  Fsin0)l  , 

me 

= - —  [— U1  (—Bo  sin  koz  —  E  cos  0)  —  U 2  (Bq  cos  koz  +  E  sin  0)]  . 


Making  the  substitutions 


u 1  =  c/3\ 

dt 

dT  =  — , 

7 


(1.6) 


(1.7a) 

(1.7b) 

(1.7c) 

(1.7d) 

(1.8) 

(1.9) 


'In  a  helical  undulator,  consecutive  magnetic  fields  are  rotated  90°  from  each  other  causing  the  electron 
to  follow  a  corkscrew  path.  The  derivation  for  electron  motion  in  a  linear  undulator  is  more  complicated. 
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where  7  =  (1  —  (3  •  (3)  1/2,  we  can  rewrite  (1.7)  as 
d'y  c 

—  = - E  (f3x  cos  ip  —  (3y  sin  ip)  ,  (1.10a) 

at  me 

77(7  At)  =  [Ecos'ip  -  f3z  (B0sin  k0z  +  Ecos'ijj)} ,  (1.10b) 

at  me 

d  c 

—  (7 (3y)  = - [-E  sin  i/j  — /3Z  (—Bq  cos  koz  —  E  sin  ip)\ ,  (1.10c) 

az  vne 

d  c 

yiilPz)  = - -[E  (f3x  cos  ip  -  (3y  smijj)  +  B0(f3x  sin  k0z  -  Py  cos  k0z)\ .  (l.lOd) 

clz  vne 

The  magnetic  held  does  no  work,  and  indeed  (1.10a)  states  that  the  electron’s  energy  is 
affected  only  by  the  coupling  of  the  radiation  electric  held  to  the  electron’s  transverse  motion. 
It  is  convenient  to  rewrite  this  equation  as 

7  =  '  E.  (1.11) 

me 

By  defining  a  transverse  velocity  (3±  =  /3xit  +  /3y y  =  (/3X,  /3y,  0),  we  may  combine  (1.10b)  and 
(1.10c)  into  a  single  expression 


dh/3±)  £ 

—— - = - \E  (1  —  /3Z)  (cos  ip,  —  sin^,  0)  +  /3*Bo  (—  sin  koz,  cos  koz,  1 

at  me 

which  for  highly  relativistic  electrons  (/3Z  — >  1)  can  be  approximated  by 


d{iP  1 

dt 


~ - [8zBq  (—  sin  koz,  cos  koz,  0)1 . 

me 


Equation  (1.13)  is  readily  integrated: 

J  rf(7 ’P±)  = 

iPx  = 

P±  = 


eBi 


0 


me 

eB0 

me2 


eBn 


ko'jmc2 

eBo\o 


/3Z(—  sin  koz,  cos  koz,  0 )dt, 
-^(—  sin  k0z ,  cos  k0z,  0)0, 
cos  koz,  sin  k^z,  0)  +  C 
(cos  k0z,  sin  k0z,  0)  +  C. 


2'Kymc2 

We  eliminate  the  constant  of  integration  by  assuming  perfect  injection,  so  that 

w°>  =  “S' 

Using  (1.2)  and  the  fact  that  for  a  helical  undulator  Brms  =  B0,  we  can  write 

K 

f3_ l  = - (cos  koz,  sin  k^z,  0) 

7 


(1.12) 

(1.13) 

(1.14) 

(1.15) 

(1.16) 

(1.17) 

(1.18) 
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or 


Px  = - cos  k0z, 

*  - 


(1.19) 

(1.20) 


Substituting  (1.18)  into  (1.11),  we  have 

eK  E 

7  =  - (cos  k0z,  sinfc0z,  0)  •  (cosi/j,  —  sin^,  0) 

7  me 

=  — —  (cos  k0z  cos  i/j  —  sin  k0z  sin  p) 

7  me 

eKE  n 

=  - cost  koz  +  ip). 

ymc 

Meanwhile,  substituting  (1.19)  and  (1.20)  into  (l.lOcl)  gives  us 


—  (7 Pz)  = - E[ - cos  k0z  cos  p  H - sinfcozsin^ 


(  K  K 

+B 0 - cos  koz  sin  k^z  H - sin  koz  cos  k^z 

V  7  7 

eK 

PPz  +  7 Pz  =  - [ E  (cos  koz  cos  i/j  —  sin  k^z  sin  ip) 

ymc 

—B0  (cos  k0z  sin  k0z  —  sin  k0z  cos  k0z)\ , 

eKE 

iPz  =  - COS (k0z  +  ip)  -  7 pz. 

7  me 


Substituting  (1.21)  into  (1.24),  we  have 


eKE  eKE 

7 Pz  =  - cos  (k0z  +  ip) - cos(k0z  +  ip)Pz, 

7  me  ymc 

eKE 

iPz  =  (1  -  Pz) - COs(k0Z +  1p), 

ymc 

eKE 

pz  =  (1  -  Pz)—- — cos{kQz  +  ip). 

7  Amc 


(1.21) 


(1.22) 


(1.23) 

(1.24) 


(1.25) 

(1.26) 


(1.27) 


It  is  useful  at  this  point  to  collect  the  electron  terms  in  k^z  +  ip  by  identifying  the  electron 
phase 

(  =  (k  +  k0)z  —  cut,  (1-28) 

referenced  to  the  combined  optical  and  undulator  fields.  This  leads  to 


k0z  +  ip  =  (  +  p. 


(1.29) 


(  =  (k  +  k„)y-u 

at 

=  (k  +  k0)pzc  -  u. 


(1.30) 
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Solving  for  /3Z  gives 


P,  = 

Pz  = 


(  +  u 

(k  +  k0)c 

c 


(k  +  k0)c 

The  expression  for  f3z  may  be  simplified  as  follows:  for  laser  light  we  require 

27 rc 


u  = 


A 


Meanwhile,  the  change  in  electron  phase  over  the  course  of  the  undulator8  is 


•  27 rc  9  _  j 

C  ~  ~  10  s 

A0 


In  other  words,  since  typically  Ao  3>  A,  then  oj  (  and  (1.31)  becomes 

0  =  —^ _ . 

(k  +  k0)c 

Using  (1.29),  (1.32),  and  (1.35),  we  may  now  rewrite  (1.27)  as 

c  eKE 


(1.31) 

(1.32) 

(1.33) 

(1.34) 

(1.35) 


UJ 


cos(C  +  0), 


(k  +  ko)c  7 2,mc  [  (k  +  ko)c\ 

eKE 

C  =  -7—;  [(k  +  k0)c  -  u]  cos(C  +  0) 


7  *mc 
eKEk0 
7  2m 


cos  (C  +  0). 


(1.36) 


(1.37) 


We  can  normalize  the  evolution  of  the  electron  motion  to  the  length  of  the  undulator  by 

o 

defining  a  dimensionless  time  t  =  ct/L.  We  adopt  the  notation  (=  dQ / dr  and  we  write 


c  = 


eKEkpL 2 
7  2mc2 
2tt  eKEL2 


cos(C  +  </>) 
cos(C  +  4>) 


72mc2A0 
27t  NeKEL 

- - - - - COs(C  +  0), 

7 


(1.38) 


where  TV  =  L/Ao  is  the  number  of  undulator  periods.  Defining  the  dimensionless  optical 
field 

_  2tt NeKEL 

«  = - 5 - o — ’  (1-39) 

7  zmcz 


8To  achieve  gain,  the  electrons  must  “bunch”  in  phase  as  they  traverse  the  undulator.  This  is  discussed 
further  in  the  next  two  sections. 
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we  have  the  pendulum  equation 


C=  l“l  cos«  +  <t> ) 


(1.40) 


Figures  1.4  and  1.5  show  the  results  of  a  simple  simulation  of  an  operating  FEL.  The  electron 
trajectories  in  phase  space  are  indeed  of  a  pendular  nature:  the  paths  are  elliptical  (closed 
paths,  as  in  Figure  1.4  -  the  pendulum  swings  back  and  forth),  and  there  exists  a  separatrix, 
a  boundary  outside  of  which  change  in  electron  phase  (  does  not  reverse  sign  (open  paths, 
as  in  Figure  1.5  -  the  pendulum  goes  “over  the  top”).  These  figures  will  be  discussed  more 
in  the  following  sections. 


Figure  1.4:  Electrons  moving  in  phase  space,  within  the  separatrix  (u0  =  0). 
Although  there  is  some  evolution  of  the  optical  phase,  there  is  no  significant 
gain  on  resonance.  (The  shading  of  electrons  darkens  with  increasing  r.) 


Figure  1.5:  Electrons  moving  outside  the  separatrix  (no  =  2.6).  There  is  gain 
due  to  favorable  bunching  of  the  electrons  in  phase  space. 


2.  Evolution  of  the  Optical  Mode 


The  fields  generated  by  the  oscillating  electrons  must  satisfy  the  inhomogeneous  Maxwell 
equations 


da^afS 


c 


(1.41) 


where 


^af} 


(  0  —  Ecosip  E  sin  ip 

E  cos  ip  0  0 

— Esin'ip  0  0 

\  0  —  E  cos  ip  Esin'ip 


°  \ 

E  cos  ip 

—E  sin  ip 

0  / 


(1.42) 


is  the  radiation  part  of  (1.6),  and  J“  =  (cp,  J)  is  the  electron  beam  4-current.  We  seek  a 
4-vector  potential  Aa  =  ($,  A)  satisfying 


=  iTA0  -  a13 A", 


(1.43) 


and  the  Lorenz  condition  daAa  =  0,  for  which  (1.41)  becomes  the  wave  equation 

da  (daA<3  -  d^Aa)  =  —  J13, 
v  '  c 

A  nr 

dadaAp =  — J P, 

c 

□A"  =  —  Ja, 
c 

where  O  =  dada  is  the  d’Alembertian  operator.  A  suitable  choice9  is 


(1.44) 


which  gives 


Aa  = 


(  0  \ 


E  .  , 

—  sm 
k 

E  , 

—  cos  ip 
k 


V  0  J 


daAa 


d  E  .  8  E  ,  n 

—  —  sm  ip  +  —  —  cos  ip  =  0 
ox  k  ay  k 


(1.45) 


(1.46) 


9This  choice  is  suitable  because  it  corresponds  to  a  helically  polarized  wave.  A  general  expression  for  Aa 
in  terms  of  the  electron  4- velocity  is  available  in  the  form  of  the  Lienard-Wiechert  potentials  (see  pp.  661-665, 
J.  D.  Jackson,  Classical  Electrodynamics,  Wiley,  New  Jersey,  1999).  Our  choice  is  made  with  the  benefit  of 
hindsight  and  is  not  unique. 
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and,  using  ca  =  kc, 


9x4° 

<9x4° 

=  0, 

dx0 

dx0 

dA 1 

<9x4° 

E  d 

<9x0 

dx\ 

kcdt 

<9x42 

<9x4° 

E  d 

(9a;o 

dx2 

kcdt 

(1.47a) 

(1.47b) 

(1.47c) 


satisfying  (1.42).  We  now  write  (1.44)  as  two  second-order  partial  differential  equations 

1  d2  E  .  d 2  E  .  4t r  T  ^  , 

?wlsm'4’-d?lsm,p  =  TJ*’  (1'48a) 


Id 2  E  ,  d 2  E  ,  4vr  T 


For  E(z,t)  and  4>(z,t)  slowly  varying,  that  is 


d2E  ,  dE 
Ik2  <  k~dz  ’ 


d2E  dE  d2(p  d(f>  d2<f  d<j) 
W  ^  u~dt  1  dk  ^  k~d~z  ’  di2  ^  udi 


we  keep  only  the  boxed  terms  in 

1  92£  .  ,  2  dEdd> 


2  9E50  2dE 

- —  COS  ip - cos  ip 

cue  dt 2  c  dt 


E  d2<p  ,  E  (d<\> 

- —  cos  ip - - 

LUC  dt 2  LUC  V  dt 


2 E  d<p  luE 

- sin  ip - sin  ip  — 

c  dt  c 


d2E  .  dE  ,  2  cdEd<j> 

-  -7— sin  ip  -2— —cos  ip - — r-cos^ 

1  (72:  a; 


cB92B  ,  caE  ,  30  ,  uE(d<t>\2 

- —  cos  ip  -\ - sin^H-  21s — sin  ip  H - I  —  sm  ip 

lu  dz 2  c  dz  c  \dz  J 


(1.48b) 


(1.49) 


(1.50a) 


leaving 


1  a2E  ,  2  dEd<f>  ,  2  8E  ,  E  d2<t>  .  ,  £  /9<A2 

—  “H72"cos^ - ^72_sln,/’  +  "  ~a7  sm ^ - - WT  cosV> 

cdc  o;c  arz  c  ot  luc  otz  luc  \ot  J 

2 E  <90  1  wE  cd2E  dE  2 cdEd<f> 

+ - cos  ip - cos  ip - —  cos  ip  -\-  2 - sin  ip  H - —  sin  ip 

c  dt  c  lu  dz2,  dz  lu  dz 2 


cE  d2E 


sin  i\)  H - cos  ip  +  2  E  — -  cos  ip  + 

dz2  c  dz 


luE  /  d(f> 


(1.50b) 


^  ( 1  d(j>  d(j)\  (IdE  dE\  An 

2E  [cm  +  T;)  ^  ~ 2  [cm  +  J7)  ^  =  —Jx'  (L51a) 

fldcb  dcb\  f  l  dE  dE\  47t 

2£  (c  al +  af) cos^  + 2  (c  aF  +  sj) sin^  =  TJ>-  (L51b) 

Defining  Jj_  =  (JX1  Jy,  0),  this  can  be  written 

{IdE  d  E\  ,  ,  x  fldcb  dcb\  4tt  , 

2  UI+&  (_“,l’SinWI  +  2i!U*+fc  <Si11  A  “s  7-.  o)  =  -Jx  .  (1.52) 


and  identifying  the  orthogonal  unit  polarization  vectors 


we  can  write 


The  operator 


Ci  =  (—  cos-0,  sin-0,  0) 
e2  =  (sin?/?,  cos-0,  0) 


d  E  d E\  . 

9f  +  caF  =  2lrJ±'ei' 


£(s+7j)  =  2irJl"2- 


d  d  d  dz  d 

dt  +  C  dz  dt  dt  dz 

d 


(1.53a) 

(1.53b) 


(1.54a) 


(1.54b) 


(1.55) 


is  the  rate  of  change  in  coordinates  traveling  with  the  evolving  light  wave,  and  we  write 


dE  2ttL 

~r~  =  - J-L  ‘  ei ) 

dr  c 

d(f)  _  2nL 

E —  —  - 

dr  c 

The  transverse  current  is  the  sum  of  individual  electron  currents,  that  is 


(1.56a) 

(1.56b) 


J±  =  -ec y^j3±  0 ! x  -  r4), 


(1.57) 


or,  using  (1.18), 


<5(x  —  r  0  (cos  k0z,  sin  k0z,  0) . 


(1.58) 


Equations  (1.56)  are  then 


dE  2nKeL 


y]  <5(x  —  rj)(—  cos  k0z  cos?/?  +  sin  k0z  sin-0),  (1.59a) 


U  2vr  KeL^., 

E—  =  -  >  o(x  —  rj (cos  Ko^rsin -0  +  sin  k^z  cos -0), 

(i.T  'V  *  ^ 


(1.59b) 


which  may  be  written 


2n  KeLnP 


dr  7  Jo 

d(j>  2n KeLne  f2n 

dr  7  J o 


COs(C  +  0)  — , 
sin(C  +  0)^, 

Z7T 


(1.60a) 


(1.60b) 
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where  ne  is  the  electron  number  density,  and  individual  electrons  are  identified  by  their 
initial  coordinates  in  phase  space  (Co,  z'o)-  All  electrons  are  assumed  to  have  the  same  initial 
energy  vq.  These  may  be  combined  into 


d,E  d(f> 

ih  + 


27 tK  eLne 

7 

2irKeLne 

7 


i  sin(C  +  (ft) 


d(o 

2vr’ 


Multiplying  both  sides  by  gives 


d  Eci4>  =  2?r  KeLne  f2w  g_ifdCo 

dr  7  J0  2vr ' 


A  compact  form  of  this  expression  is  found  by  identifying 


(1.61) 

(1.62) 


(1.63) 


a 

J 


a  e 


8n2K2e2L2Nne 
73 me 2 


(1.64) 

(1.65) 


where  |a|  is  given  by  (1.39).  Letting  an  overbar  indicate  averaging  over  all  initial  electron 
phases,  we  have 

(1.66) 


It  is  clear  that  there  is  no  change  in  optical  field  (da/ dr  =  0)  for  j  =  0.  For  gain,  the 
integral  further  requires  that  the  initial  phase  distribution  of  electrons  average  to  less  than 
zero.  Populations  of  electrons  for  which  e~ ^  ^  0  are  said  to  be  bunched  in  phase,  and 
the  manipulation  of  electron  phase  distributions  to  achieve  favorable  bunching  is  a  critical 
concept  in  FEL  design. 


3.  Interaction  and  Gain 

We  now  have  expressions  for  the  evolution  of  the  electrons  in  the  undulator  (1.40)  and 
the  light  they  create  (1.66).  To  show  how  this  light  is  amplified,  we  introduce  a  solution  to 
(1.40)  for  weak  fields  (|o|  <C  7r)  and  low  current  ( j  <  7 r).  It  is  an  expansion  in  ao,  namely10 

C(t)  =  Co  +  vqt  -  — [cos (Co  +  vqt)  -  cos(Co)  +  ^oTsin(Co)]  H - ,  (1.67) 

uo 

o 

where  v  =  C,  ao  is  the  initial  optical  field,  and  the  initial  phase  is  taken  to  be  zero.  Thus11 
Ht)  =  uo  +  —  [sin(Co  +  "or)  -  sin(Co)] 
a2  (  1 

H — Is-  -[cos(2Co  +  2 u0r)  -  cos2Co]  (1.68) 

h  l  4 

+  cos  uq t  -  1  —  uqt  sin  Co  cos(Co  +  vqt)  \  H - 


10W.  B.  Colson,  “Free  Electron  Laser  Theory”,  doctoral  dissertation,  Stanford  University,  1977. 
nIbid. 
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We  can  use  (1.67)  and  (1.68)  to  plot  the  phase  space  portrait  of  populations  of  electrons  over 
their  travel  time  through  the  undulator,  but  a  more  accurate  picture  is  achieved  numerically 
using  a  computer  simulation,  such  as  the  one  used  to  generate  Figures  1.4  and  1.5. 

The  amplification  of  the  optical  mode  is  proportional  to  the  energy  loss  of  the  electron 
beam.  Leaving  (1.68)  behind  for  the  moment,  we  can  examine  the  change  in  individual 
electron  phase  velocity  v  by  differentiating  (1.28): 


dc 

dr 


d 

v  =  -7^[(k  +  k0)z  -  ut] 


L 

c 


/  r  ,  \dz 

(k  +  k0)—  -u 
at 


=  L[(k  +  ko)(3z  —  k\, 


(1.69) 


which,  assuming  k  /c0,  may  be  approximated 


v  =  Lk(f3z  —  1),  (1.70) 

that  is 

Av  =  LkA(3z.  (1.71) 

This  is  not  a  surprising  expression:  change  in  electron  phase  velocity  is  certainly  propor¬ 
tional  to  change  in  its  longitudinal  velocity,  as  the  Lorentz  force  felt  by  the  particle  as  it 
passes  over  each  magnet  in  the  undulator  depends  on  this  velocity.  Moreover,  proximity  to 
the  resonance  condition — most  certainly  a  function  of  longitudinal  electron  velocity! — will 
ultimately  determine  the  gain.  Near  resonance,  we  may  write  (see  Equation  1.1): 

47T72 

’  =  A0(l  +  A'2)  ’ 

thus 

AnL^2Af3z 
=  vTTa7 

AnN'y2  A(3Z 
1  +  K 2 

To  put  this  expression  in  terms  of  electron  energy,  we  must  re-examine  the  relationship 
between  /3Z  and  the  relativistic  7.  In  our  helical  undulator,  where  (1.18)  implies 


(1.72) 


(1.73) 


the  definition 


P. 


2 

x 


1 

a/1  -  /3  •  /3 


(1.74) 


(1.75) 
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becomes 


1 


7  = 


'1-/31-  (32 


(1.76) 


which  may  be  written 


/  K2 

'l  ~--Pl 

7 


1  +  A'2 


(1.77) 


1  -(3i  =  — 
7 


and  so,  taking  the  differential  of  each  side, 

(3ZA/3Z  = 


A7(l  +  K 2) 


(1.78) 


(1.79) 


A&  = 


A7(l  +  A'2) 


Plugging  this  into  (1.73)  we  have 


A  u  = 


AnNA'y 


(1.80) 


(1.81) 


At  this  point  we  may  return  to  our  approximation  (3Z  — >  1  for  relativistic  electrons,  and  we 
write 


Av  =  AnN—. 


(1.82) 


For  an  electron  of  energy  Ee  =  'jmc2,  the  energy  change  as  a  function  of  the  change  in 
electron  phase  velocity  is  then 


AAe  = 

AnN 


(1.83) 


Now,  the  number  of  electrons  in  a  volume  14  of  electron  beam  is  ne14-  The  change  in  electron 
beam  energy  A E^  corresponding  to  14  is  thus 


.  „  neymc  - —  r 

=  “toir A"14’ 


(1.84) 


where  we  have  averaged  the  change  in  phase  velocity  for  individual  electrons.  Meanwhile, 
the  optical  mode  energy  E0  in  a  volume  14  is  simply 


E2 

E0  =  —VOJ 

An 

where  E  is  the  optical  field  strength.  We  now  write  the  optical  field  gain  as 

G  = 


(1.85) 


(1.86) 


nejmc2  Vb— 

~Nmr^v 

nejmc2  F  -r— 


(1.87) 
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where  the  filling  factor  F  specifies  that  fraction  of  the  optical  mode  coincident  with  the 
electron  beam.12  Using  our  definitions  for  the  optical  amplitude  (1.39)  and  beam  current 
(1.65),  this  can  be  simplified  to 

G  =  -2^-  ~Av.  (1.88) 

ao 

To  find  an  expression  for  A v,  we  return  to  our  weak  held  expansion  of  v(r)  (1.68)  and 
integrate  over  all  possible  initial  electron  phases  Co: 

Av  =  v  —  vq 


f*2n 


2n  L  ("-"°)(iCo 


(  Oo 
U) 


r»27r 


sin  (Co  +  "or)  ~  sin(Co) 


2vr 


+^f  {  ~  |  [  cos(2^°  +  2z/°r)  _  cos  2^°] 


+  cos  v0t  -  1  -  U0T  sin  Co  cos(Co  +  Vqt)  \  H - 


>dCo-  (1-89) 


Many  of  the  terms  average  to  zero,  and  so,  to  second-order  in  a0,  we  have 

2  cos  Vq t  —  2  +  u0t  sin(z/0r) 


A  v  ~  o  (2 


2z/3 


(1.90) 


Plugging  this  into  (1.88)  and  dividing  through  by  jF,  we  have 


G 

2  cos  uqt  —  2  +  vo  t  sin(i/0T) 

Jf  = 

V3 
u 0 

(1.91) 


which  describes  the  FEL  gain  along  the  undulator.  The  final  FEL  gain  at  r  —  1  (the  end  of 
the  undulator)  is  shown  in  Figure  1.6.  The  salient  feature  of  this  figure  is  its  antisymmetry, 
in  particular  its  behavior  close  to  is0  =  0,13  which  indicates  that  energy  may  be  either 
transferred  from  the  electron  beam  to  the  optical  mode  (0  <  vq)  or  vice-versa  (vq  <  0). 
Earlier  (see  p.  7),  we  remarked  that  electrons  must  bunch  in  phase  space  to  transfer  energy 
to  the  optical  mode.  For  a  functioning,  useful  FEL,  (1.91)  narrows  that  requirement  to  a 

12For  electron  beam  and  optical  mode  volumes  both  spanning  a  longitudinal  length  l,  the  filling  factor  is 
defined  as  a  ratio  of  cross  sectional  areas: 

Vb  _  7r rlX  _  Ab 

V0  7T  r^X  A0 

13It  is  interesting  to  examine  why  (1.91)  does  not  blow  up  at  vq  =  0.  After  all,  the  presence  of  in  the 
denominator  is  certainly  suggestive  of  singular  behavior  at  the  origin.  The  answer  is  found  by  expanding 
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Figure  1.6:  The  FEL  weak  field  gain  curve  plotted  at  the  end  of  the  undulator 
(t  =  1).  Energy  is  either  gained  by  the  optical  mode  or  lost,  depending  on 
initial  phase  velocity  Uq. 


population  of  electrons  possessing  an  average  phase  velocity  just  to  the  positive  side  of  the 
G/jF  axis  in  Figure  1.6. 

To  illustrate,  we  refer  once  again  to  Figures  1.4  and  1.5.  These  show  the  results  of  a 
simple  simulation  of  FEL  gain  for  an  electron  beam  with  z/o  =  0  and  uq  >  0.  In  the  first  case 
(Figure  1.4),  nearly  the  same  amount  of  energy  gained  by  the  optical  mode  from  electrons 
with  decreasing  phase  velocity  is  simultaneously  lost  to  electrons  whose  phase  velocity  is 
increasing.  In  the  second  case  (Figure  1.5),  where  u0  >  0,  electrons  bunch  in  the  desired 
region  of  phase  space:  they  are  predominantly  losing  energy  to  the  optical  mode. 


the  numerator  about  vq  =  0,  where 

2  cos  vqt  —  2  +  t'o  t  sin(r'or) 


2  -  u0t 


,,4^4 


2  2  ,  "0  1 


,4^4 
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-  2  +  V0T  (  V0T  - 


.3—3 


12 


In  the  vicinity  of  the  origin,  then,  at  r  =  0  (1.91)  becomes 

G  1 

=  - Un, 

jF  12 


confirmed  by  the  linear  nature  of  the  vq  «  0  region  of  Figure  1.6. 
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The  introduction  of  the  gain  equation  completes  our  preliminary  discussion  of  free  elec¬ 
tron  lasers.  For  convenience,  we  review  the  important  conclusions  of  this  chapter:  the  two 
principle  equations  of  basic  FEL  theory,  and  an  expression  (derivable  from  the  first)  for  gain 
in  weak  fields. 


Pendulum  Equation 

oo 

c=  |a|  cos(C  +  <t>) 

Wave  Equation 
a  =  —j 

Gain  Equation 

(weak  fields) 

G  2  cos  u0t  —  2  +  u0t  sin  (z/0t) 

jF  ul 


(1.40) 


(1.66) 


(1.91) 
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Chapter  II. 


Fundamental  Mode  Optics 


The  free  electron  laser  is  unique  among  high  power  lasers  in  that  it  operates  primarily 
in  the  fundamental  mode.  In  contrast,  high-power  solid  state  and  chemical  lasers  deliver  a 
beam  comprised  of  many  higher-order  mode  (HOM)  components.  Over  a  given  propagation 
range,  these  HOMs  lead  to  more  diffraction  loss  than,  for  example,  a  purely  Gaussian  beam. 
While  such  a  beam  characterizes  the  ideal  laser,  we  will  show  in  the  next  chapter  that 
beams  obtained  from  simulations  of  FELs  appear  to  contain  a  few  HOMs,  and  it  is  our  goal 
to  construct  a  diagnostic  tool  to  evaluate  the  mode  content  of  a  given  laser  beam  whose 
amplitude  and  phase  are  known  at  every  point. 

Such  comprehensive  knowledge  of  beam  makeup  is  possible  in  the  laboratory  with  much 
effort,  but  a  simulation  provides  all  such  information  a  priori  and  allows  a  comprehensive 
understanding  of  the  system,  and  the  ability  to  alter  that  system,  without  the  high  cost 
of  a  laboratory  experiment.  These  types  of  simulations  have  been  benchmarked  against 
actual  experiment  and  have  accurately  predicted  experimental  results.  Modal  analysis,  on 
the  other  hand,  has  not  been  benchmarked,  and  this  analysis  provides  a  starting  point  for 
such  a  comparison. 

A.  Propagation  -  Analytic  Description 

Although  we  intend  to  approach  the  FEL  numerically,  we  are  interested  in  an  analytic 
solution  to  the  wave  equation  in  order  to  understand  the  nature  of  an  evolving  wavefront 
inside  a  laser  cavity.  For  the  moment,  we  consider  the  optical  mode  independently  of  its 
interaction  with  the  electron  beam. 

1.  The  Parabolic  Wave  Equation 

The  homogeneous  part  of  (1.44)  can  be  written 

(v2-lA)A(x,t)=0.  (2.1) 

Expressions  solving  this  equation  describe  waves  propagating  in  free  space  (or  other  homo¬ 
geneous,  sourceless  media).  Using  the  same  vector  potential  as  in  the  previous  chapter,  we 
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have 

=  ^  (sin-0,  cos0,  0),  (2.2) 

k 

where  e  =  (— i,  1,  0)  is  the  polarization  vector,  a  =  kz  —  out  is  the  carrier  wave  phase,  and 
£(x,  t)  =  E0x,  t)el<^x’^  is  comprised  of  an  amplitude  E  and  phase  0  now  with  an  explicit 
transverse  dependence  (remember  that  before  we  used  E  =  E(z,t)  and  0  =  <f>(z,t );  see 
page  10).  To  begin  substituting  this  into  the  wave  equation  (2.1),  we  take  two  derivatives  in 
z  and  t: 


A(x,  t)  =  3? 


£(x, 


dA 

dz 

o2a 

dz2 

dA 

~dt 

82A 

~W 


^KI+!K)ei“e}' 

kS+2tk^-k2£]e'^’ 

3?  <!  -(^  -  2iu^-  -  u2S  )  e*Qe 

dt 


k  \  dt2 


(2.3) 

(2.4) 

(2.5) 

(2.6) 


As  before,  for  slowly- varying  Efx.,t)  and  0(x,  t),  variations  in  amplitude  and  phase  arising 
from  diffraction  in  the  direction  of  propagation  are  small  compared  to  those  of  the  wave 
itself,  that  is 


d2E 

,  dE 

dz2 

•C 

k — 

dz 

d2E 

dE 

dt 2 

•C 

uj— — 

dt 

d2<f> 


dz2 


•C 


k9^ 

kdz 


d2cj) 
dt 2 


•C 


(2.7) 


These  negligible  terms  comprise  the  second  longitudinal  space  and  time  derivatives  of  £, 
which  we  eliminate  to  obtain 


d2£  d2£ 


dx2  dy 2 


,  dS  ,  9  A  1  (  d£  2  r> 
2 ik  — - k2£ - -  I  —2 iuj— - uj2£ 


dz 


dt 


eiae  =  0. 


(2.8) 


This  simplifies  to 


V2£  +  2ik  [^-  +  ~ 

dz  c  dt 


eiae  =  0, 


where  we  have  used  u  =  kc  and  introduced  the  transverse  Laplacian 

d2  d2 


= 


+ 


dx 2  dy2 


Multiplying  by  cue  *a,  the  scalar  equation  is 

2  „  ,  (d£  d£  . 

cV±£  +  2 ik  ^—+c—)-0, 

and  we  may  once  again  use  (1.55)  to  write  this  as 


1 _ r,  ik  d  .  „ 

2V±  +  J(h  =  0- 


(2.9) 


(2.10) 


(2.11) 


(2.12) 
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With  propagation  now  expressed  in  terms  of  r,  we  would  like  to  make  the  diffraction  term 
dimensionless  also.  The  1/e  radius  w o  of  the  initial,  fundamental  mode  of  the  laser  is  given 


by 


z0X 

71 


(2.13) 


where  zq  is  the  Rayleigh  length,  or  the  characteristic  distance  over  which  the  beam  cross 
section  increases.  Defining  the  dimensionless  transverse  coordinates  as 


so  that 


81 2  _  k  d2  d2 
dx2  2  L  dx2  ’  dy 2 


and  adopting  the  dimensionless  transverse  Laplacian 


k  d2 
2Ldtp’ 


d2  d2 
dx 2  dy 2 1 


we  now  write  the  dimensionless  parabolic1  wave  equation 


(2.14) 

(2.15) 

(2.16) 

(2.17) 


2.  Solving  the  Dimensionless  Parabolic  Wave  Equation 

To  solve  (2.17),  we  attempt  a  trial  solution  of  the  Gaussian  form2 

S(x,y,r)  =  £op(r)  exp 


(  x2  +  y2\ 

V  ?(r)  )’ 


(2.18) 


where  p(r )  and  q(r )  are  as  yet  undetermined  complex  functions  describing  the  evolution  of 
amplitude  and  width  over  dimensionless  time.  We  further  impose  the  initial  condition 


£(x,y,0)  =  £0  exp 


^  x2  +  y2^ 


(2.19) 


1The  terms  paraxial  and  parabolic  may  both  be  used  to  describe  (2.17).  It  is  paraxial  in  that  it  accurately 
describes  wave  behavior  close  to  the  axis  of  propagation.  It  is  parabolic  in  that  it  fits  into  a  class  of  partial 
differential  equations  of  the  general  form 


d2u  d2u  „d2w  du 

Adx?  +  2Bdxdy  +  +  Ddx  +  Edy 


+  Fu  =  R(x,  y), 


where 

=  AC  —  B2  =  0. 

This  class  of  equations  includes  the  Schrodinger  equation,  with  which  the  parabolic  wave  equation  shares 
many  properties.  In  contrast,  the  original  wave  equation  is  hyperbolic,  that  is  AC  —  B 2  <  0. 

2 A  spherical  wave  is  a  known  solution  to  the  wave  equation.  A  Gaussian  wavefront  may  be  described  as 
a  paraxial  spherical  wave,  an  approximate  solution  to  the  wave  equation.  It  is  also,  in  fact,  an  exact  solution 
to  the  paraxial  wave  equation. 


A  B 
B  C 


21 


where  z0  =  Wq  is  the  dimensionless  Rayleigh  length  and  the  central  held  amplitude  is  £0  = 
£(0,0,0).  Substituting  (2.18)  into  (2.17),  we  take  the  appropriate  derivatives  in  x,  y ,  and  r 
(letting  r2  =  x2  +  y2  and  suppressing  the  argument  r): 


d2£ 
dx2 
d2£ 
dy 2 


£0pe-f2/q, 

(U 

-)  £0 pe-f2'\ 

q) 


(2.20) 

(2.21) 


d£  (l  dp  r2dq\  _r-2/ 
=  {pTr  +  ^Tr)£oPe  ’ 


(2.22) 


and  write 

(4--) 

\  r  qj 


£0  pe-f2/q  + 


tU^  +  MilfT  +  ^r)saPe 


D-f2/q 


0.  (2.23) 


which  may  be  written3 


For  this  equation  to  hold  for  all  r  (that  is, 
must  be  identically  zero,  thus 

dq 

dr 

dp 

dr 

It  is  easy  enough  to  solve  (2.26)  for  q(r): 

q(r) 


(2.24) 

(2.25) 

the  expressions  in  parentheses 

(2.26) 

(2.27) 

go  +  iT-  (2.28) 


Dividing  through  by  4 £ope  r2^q  and  collecting  terms,  we  have 

7,2  1  .  fldp  r 2  dq ' 


q-  q 


_  T  +  M  “73T  +  ~^~7Z  ~ 


p dr  q2  dr 


-d/  +  -)?  =  o- 

pdr  q) 
for  all  x  and  y), 


= 


We  are  not  surprised  to  hnd  that  it  is  complex,  as  it  must  describe  not  only  the  evolution  of 
the  transverse  extent  of  the  wave,  but  also  the  evolution  of  the  phase,  which  must  oscillate. 
From  our  initial  conditions,  we  hnd 

go  =  zo,  (2.29) 

also  anticipated.  Using  this  result  and  substituting  (2.28)  into  (2.27),  we  have, 


dp  ip 

dr  So  +  it 

P 

iz0  -  t’ 

3To  get  from  (2.24)  to  (2.25),  multiply  through  by  —  iq2,  then  collect  terms  in  f2  and  q2 . 


(2.30) 


22 


which  we  solve  for  p(r): 


I o  iz0  -  r 


=  ln(iz0  -  t') 


=  ln(iz0)  —  In  (izQ  —  r) 


=  In 


1  +  i~z~  I  ’ 

Zq) 


(2.31) 


p(t)  = - r- 

1  +i— 
z0 


(2.32) 


What  does  it  mean  for  p(r)  to  be  complex?  Our  solution  to  the  parabolic  wave  equation  is 


x2  +  y2 


S{x,y,r)  = - f  exp  -  ,  .  (2.33) 

1  +  i—  V  zo  +  it  J 

Z0 

Once  again  letting  r2  =  x2  +  y2,  we  elucidate  this  expression  by  making  the  following 
definition: 

3  =  1  +  ip,  (2.34) 


where 


and  we  write 


F(-  ~  v  So  (  f 2 
£{x,y,T)  =  —  exp  -  — 

3  V  2o3 


£o  f  r2  3* 

—  exp  - - - 

3  V  zo  33 


So  F  n  ■  )  r~2 

=  -exp  -(l-tp)  — 
3  L  z o33 


,p2  ~2 

- - h  ip- - 

2o33*  -033* 


£o  (  f 2  \  (■  f2 

—  exp  -- — -  exp  ip- — - 

3  V  ^o33*  /  V  zo33 


Now,  the  complex  number  3  has  modulus 


(2.34) 


(2.35) 


(2.36) 


I3I  —  \/l  +  P2 ) 


(2.37) 
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and  argument 
so  we  may  write 

and  (2.36)  becomes 
S(x,y,r)  = 


arg3  =  arctan  //, 

j  g— i  arctan  fi 


3  y/l  +  yU2  ’ 


£n 


:  exp 


\/l  +  /P 

Defining  the  dimensionless  beam  width 


W^W)) exp  U  i^tt^ 


.  —i  arctan  /  /. 


W 


—  \/ A)33*  —  x/MTTT)  —  \  Zq  (  1  +  ^  ) , 


we  arrive  at 


£ (5, 1/,  t)  =  £0 -^r  exp  (  -  r  J  exp  (  i  _/  T  9  J  exp  (  -i  arctan  P  )  , 

w(r)  V  w2{r)J  \z2  +  t2J  V  -o/ 


where 


w0  =  w(Q)  =  \f%o. 

This  is  the  solution  for  the  dimensionless  parabolic  wave  equation.  Using 


0(r,  r)  = 


r2r 

zl  +  r2 


r 

—  arctan  — , 

z0 


we  write  finally 


_  C  W°  „~r2 /w2(t) ci<t>(f,T) 


£(x,y,r)  =  £0—e 


(2.38) 

(2.39) 

(2.40) 

(2.41) 

(2.42) 

(2.43) 

(2.44) 

(2.45) 


B.  Properties  of  the  Gaussian  Beam 

We  have  identified  a  physical  basis  for  each  of  the  factors  in  (2.45);  these  are  reviewed 
in  Table  2.1.  Figure  2.1  illustrates  a  beam  diffracting  from  a  waist  at  r  =  0  in  the  ±r 
directions,  identifying  the  important  parameters  used  in  the  Gaussian  beam  expression.  We 
can  study  the  behavior  of  the  beam  described  by  this  expression  by  looking  at  its  variables 
in  their  limits.  As  we  have  required,  at  r  =  0, 

S(x,y,0)=Soe-f2/*°.  (2.46) 

As  r  increases,  so  does  w:  the  beam  expands.  For  a  given  r,  the  overall  value  on  the  right 
hand  side  of  (2.45)  decreases  with  increasing  r:  the  field  strength  attenuates. 

Figure  2.2  shows  a  Gaussian  beam  propagating  through  free  space.  The  expansion  and 
attenuation  of  the  beam  are  evident.  Also  evident  is  that  the  Gaussian  beam  retains  its  func¬ 
tional  form:  there  is  a  scaling  of  the  overall  structure,  but  the  wavefront  remains  Gaussian 
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Table  2.1:  Symbols  Used  in  the  Dimensionless  Parabolic  Wave  Equation 


Symbol 

Definition 

Description  (dimensionless) 

X 

x{jL 

orthogonal  transverse  coordinate 

y 

orthogonal  transverse  coordinate 

r 

1  7 r  r 

rV  A L  =  ^ 

transverse  coordinate 

ct 

time 

T 

L 

W 

h  0 + f ) 

beam  width 

w  0 

w{  0)  = 

beam  width  at  beam  waist 

Zo 

Zo 

L 

Rayleigh  length 

after  propagation.  This  behavior  is  unique  to  unmixed,  individual  component  solutions  of 
the  paraxial  wave  equation,  but  as  we  will  show  in  Chapter  IV.,  the  equation  allows  solutions 
beyond  just  the  fundamental  Gaussian.  Moreover,  we  will  show  that  an  arbitrary  laser  (for 
example,  a  laser  with  a  top-hat  or  annular  or  just  very  “messy”  cross-sectional  wavefront 
structure)  experiences  the  inclusion  of  higher-order  modes — that  is,  higher-order  solutions 
to  the  wave  equation — as  a  loss  mechanism4  during  propagation. 

The  above  claim  demands  to  be  addressed,  but  to  further  motivate  our  discussion  of 
higher-order  modes,  we  first  present  the  results  of  past  research  in  the  field  of  high-power 
free  electron  lasers. 


4This  is  not  loss  in  the  sense  of  absorption  or  scattering,  but  rather  with  respect  to  the  amount  of  energy 
intended  to  reach  a  target  area. 
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Figure  2.1:  A  diagram  of  the  optical  mode  showing  the  relevant  (dimension¬ 
less)  physical  meanings  of  the  symbols  in  Table  2.1. 


Figure  2.2:  The  beam  described  in  Figure  2.1  on  the  interval  0  <  r  <  1. 
The  initial  Gaussian  wavefront  (left)  propagates  through  a  uniform  sourceless 
medium  (center),  diffracting  into  a  wider  Gaussian  of  lower  peak  amplitude 
(right). 
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Chapter  III. 


Motivation:  Free  Electron  Laser 
Simulations 


A  computer  simulation  of  laser  beam  propagation  through  a  sourceless  medium  was  used 
to  generate  Figure  2.2.  Such  a  simulation  implements  a  relatively  simple  numerical  method 
for  solving  the  parabolic  wave  equation.  More  relevant  simulations  (to  wit:  the  inclusion  of 
the  radiation  source — the  electron  beam!)  are  the  work  of  the  Free  Electron  Laser  Research 
Group.  Naval  Postgraduate  School  FEL  simulations  have  been  in  use  for  fifteen  years  and 
have  been  useful  tools  in  the  construction  of  new,  more  powerful  FELs.  Recently,  ever 
increasing  computing  speed  and  power  and  advances  in  simulation  technology  (better  codes) 
have  resulted  in  greater  fidelity,  and  more  robust  analysis  of  the  evolving  FEL  system.  The 
inclusion  of  expanding  coordinates  and  methods  of  characterizing  optical  beam  quality  are 
among  the  improvements  found  in  the  following  work,  completed  over  the  last  two  years. 

A.  Multiple  Modes 

Under  optimum  conditions,  the  free  electron  laser  operates  in  the  fundamental  mode.  In 
order  to  create  the  laser  beam,  the  gain  medium — a  beam  of  free  electrons — passes  through 
the  undulator,  and  thus  some  portion  of  the  optical  cavity.  The  complicated  boundary 
conditions  introduced  by  the  presence  of  a  region  of  swiftly  moving  charges  within  the  re¬ 
gion  of  propagation  argue  strongly  for  the  use  of  numerical  methods.  But  one  might  guess 
that  the  simple  evolution  of  an  initially  Gaussian  wavefront  like  that  in  Figure  2.2  is  made 
more  complicated  by  the  presence  of  the  electrons,  and  in  fact  simulations  that  include  the 
electron  beam  show  that  an  initially  Gaussian  FEL  beam  wavefront  can  assume  an  other- 
than-Gaussian  transverse  structure. 

1.  Multi-Mode  Simulations  of  a  Short-Rayleigh  Length  FEL 

Simulations  offer  the  ability  to  test  variations  in  free  electron  laser  design  parameters 
such  as  undulator  length,  electron  bunch  charge,  and  Rayleigh  length,  before  such  designs 
are  realized.  For  a  weapon-class  FEL,  two  main  objectives  are  maximizing  output  power  and 
minimizing  the  possibility  of  damage  to  resonator  optics.  In  2003,  the  FEL  Group  simulated 
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the  performance  of  a  proposed  1  //rn  high-power  FEL  with  a  Rayleigh  length  an  order  of 
magnitude  shorter  than  typical  extant  operating  FELs.1  The  three-dimensional  simulation 
(in  x,  y,  and  r)  used  self-consistent  FEL  pendulum  (1.40)  and  wave  (1.66)  equations  to 
follow  the  evolution  of  weak  optical  fields  to  steady  state,  allowing  for  the  development  of 
transverse  structure  in  the  wavefront.  For  simulations  of  this  type,  the  value  of  interest  is 
the  extraction  rj,  defined  as 

extracted  optical  power 
initial  electron  beam  power 


More  specifically,  there  are  two  expressions  for  rj  in  basic  FEL  theory:2  one  for  the  low  gain 
regime  (j  <  1), 

V  ~  rU  (3-2) 


and  one  for  high  gain  (j  3>  1), 


n 


2N 

Cm 

8  N  ' 


(3.3) 


Both  are  approximate,  as  is  the  definition  of  the  cut-off  between  high  and  low  gain. 

In  the  simulations,  design  parameters  are  independently  varied  to  determine  their  opti¬ 
mum  values  with  rj  as  a  gauge.  The  simulations  must  determine  the  optimum  resonance 
parameter 

uo  —  L[(k  +  k0)f3z  o  —  k\,  (3-4) 


where  L  is  the  undulator  length,  k  is  the  optical  wavenumber,  ko  is  the  undulator  wavenum¬ 
ber,  and  f3z o  is  the  initial  electron  velocity  (see  Chapter  1  for  details  of  these  parameters). 
An  actual  FEL  operates  at  the  optimum  i/0,  hireling  it  by  automatically  varying  A  =  2n/k 
until  a  fixed  value  of  u0  dominates  by  mode  competition. 


2.  Results 

Each  individual  simulation  run  results  in  output  like  that  shown  in  Figure  3.1.  A  series 
of  runs  are  inspected  to  determine  the  optimum  ty  for  a  given  parameter,  say  charge  q  per 
electron  bunch.  The  plot  of  extraction  rj  versus  q  at  each  optimum  u0  shown  in  Figure  3.2  is 
thus  the  result  of  more  than  one  hundred  simulation  runs.  Figure  3.2  demonstrates  the  steady 
increase  in  extraction  with  increasing  bunch  charge.  The  simulation  assumes  a  constant 
longitudinal  emittance 

£i  =  A74,  (3.5) 

where  Ay  is  the  rms  spread  in  the  electron  energies  and  lb  is  the  rms  electron  bunch  length. 
In  a  real  FEL  (and  in  later  simulations),  emittance  increases  with  bunch  charge,  reducing 
extraction  as  the  beam  quality  degrades.  But  the  effect  is  estimated  as  small  here.3 

1W.  B.  Colson,  et  al,  “Multi-mode  Simulations  of  a  Short-Rayleigh  Length  FEL”,  Proceedings,  25th 
International  FEL  Conference,  Tsukuba,  Japan,  September  2003. 

2W.  B.  Colson  in  W.  B.  Colson,  C.  Pellegrini,  A.  Renieri  (Eds.),  Laser  Handbook,  Vol.  6,  North  Holland, 
Amsterdam,  1990. 

3For  a  concise  description  of  experimental  FEL  performance  as  it  relates  to  particular  phenomena,  such 
as  electron  beam  emittance,  see  Chapter  3  of  B.  W.  Williams,  Jefferson  Lab  Free  Electron  Laser  10  kW 
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Figure  3.1:  Output  from  a  2003  simulation.  Each  output  figure  displays  8 
panels  showing,  clockwise  from  top  left:  (1)  a  plot  of  the  mode  evolution  over 
many  passes  through  the  undulator,  (2)  a  cross  section  of  the  laser  beam,  (3) 
a  list  of  input  and  calculated  parameters,  (4)  a  plot  of  output  beam  power 
evolution,  (5)  a  plot  of  gain  evolution,  (6)  the  final  phase  space  portrait  of  the 
electrons,  (7)  an  historical  plot  of  electron  phase  space  evolution,  and  (8)  a  plot 
of  the  mode  evolution  over  the  final  pass  through  the  undulator. 


The  results  of  the  2003  runs,  which  were  presented  at  the  International  Free  Electron  Laser 
Conference  in  Tsukuba,  Japan,  included  determinations  of  optimum  number  of  undulator 
periods  (N  =  14,  see  Figure  3.3),  optimum  normalized  beam  waist  radius  (a  =  0.12,  see 
Figure  3.4),  and  constraints  on  normalized  Rayleigh  length  (z0  <  0.06,  versus  the  typical  0.5 
in  use  at  the  time,  see  Figure  3.6)  and  mirror  transmission  losses  (>50%)  for  the  proposed 
laser. 

Figure  3.3  reflects  two  important  FEL  design  concepts.  The  first  is  that  too  few  undulator 
periods  result  in  a  non-functioning  FEL,  indicated  by  the  rapid  drop  off  in  extraction  for 
N  <  10.  This  is  evident  from  our  definition  of  the  dimensionless  current  density 
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(3.6) 


Since  gain  G  is  proportional  to  current  j,  we  ultimately  have 


G  oc  j  oc  N3. 


(3.7) 


Upgrade  -  Lessons  Learned ,  Technical  Report  NPS-PH-05-001,  Naval  Postgraduate  School,  Monterey,  May 
2005. 
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Figure  3.2:  Plots  of  FEL  performance  as  a  function  of  individual  design 
parameters,  in  this  case  bunch  charge  q ,  are  each  the  result  of  many  simulation 
runs. 


Meanwhile,  the  dependence  of  extraction  on  number  of  undulator  periods  initially  obeys 
(3.2),  that  is 

V  oc  (3.8) 

countering  the  trend  due  to  increasing  gain.  In  the  high  gain  regime,  extraction  appears  to 
have  no  dependence  on  N: 


V 


m 

8  N 


m 


Un2IOe2X^ne 

j^mc2 


I  tx2K 2e2Aone 


(3.9) 


12873mc2 

and  indeed  the  curve  levels  off  once  maximum  gain  is  achieved,  around  N  —  15  for  this  laser. 

Figure  3.4  illustrates  the  effects  of  two  important  FEL  concepts:  filling  factor  F  and 
(normalized)  transverse  emittance  en,  given  by 


£n  =  rrbdb, 


(3.10) 


where  7  is  the  electron  beam  energy,  rb  is  the  rms  beam  radius,  and  9b  is  its  rms  angular 
spread.  In  general,  low  emittance  is  desirable,  and  the  value  can  only  increase  as  the  beam 
travels  around  the  FEL  transport  lattice.  Considerable  pains  are  taken  to  suppress  emittance 
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Figure  3.3:  Plot  of  extraction  ij  as  a  function  of  number  of  undulator  periods 
N.  While  the  curve  is  relatively  flat  beyond  N  =  15,  a  guiding  principal  in 
FEL  design  is  using  just  enough  undulator  periods  to  get  a  working,  efficient 
laser. 


growth,  but  it  is  never  zero,  as  evidenced  when  the  electron  beam  is  focused  at  the  center  of 
the  undulator.  Since  £n  must  at  least  be  conserved,  the  narrower  the  waist  of  the  electron 
beam,  the  greater  the  angular  spread.  Too  narrow  a  waist,  and  some  electrons  fall  outside 
the  optical  mode  at  the  ends  of  the  undulator  due  to  a  large  angular  spread,  as  in  the  top  of 
Figure  3.5.  On  the  other  hand,  the  electrons  in  a  wider  beam  interact  with  a  lower  average 
optical  field  in  the  middle  of  the  undulator,  as  in  the  bottom  of  the  figure.  Increasing  F  by 
increasing  Af,  may  result  in  some  electrons  falling  outside  the  optical  mode,  as  illustrated  in 
Figure  3.5.  These  electron  still  wiggle,  but  the  radiation  they  produce  is  not  reinforced  by 
the  optical  resonator. 

The  important  result  from  Figure  3.6  is  that  for  a  small  increase  in  extraction  achieved 
by  increasing  the  Rayleigh  length  zq,  there  is  a  considerable  increase  in  optical  intensity  on 
the  resonator  mirrors.  This  is  of  particular  concern  at  present.  Since  vulnerability  of  the 
mirror  limits  the  allowable  power,  material  limitations  make  the  value  of  mirror  intensity  a 
limiting  factor  in  the  development  of  weapon-class  free  electron  lasers. 
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Figure  3.4:  Plot  of  extraction  7/  and  optimum  initial  phase  velocity  uq  as 
a  function  of  normalized  electron  beam  waist  radius  a  for  fixed  transverse 
emittance  en.  There  is  an  optimum  overlap  between  the  electron  beam  and  the 
optical  mode. 


I a(0,  y,x ) I 


Y 


Figure  3.5:  The  overlap  of  the  electron  beam  and  the  optical  mode.  A  nar¬ 
rower  electron  beam  means  larger  emittance  sn  and  a  lower  number  of  inter¬ 
acting  electrons,  yet  a  wider  beam  creates  an  interaction  region  with  a  lower 
average  optical  field. 
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Figure  3.6:  Plot  of  extraction  7/  and  mirror  intensity  as  a  function  of  nor¬ 
malized  Rayleigh  length  zq.  For  longer  Rayleigh  lengths,  a  small  increase  in 
extraction  is  accompanied  by  a  large  increase  in  (potentially  harmful)  optical 
intensity  on  the  resonator  mirrors. 
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B.  Optical  Mode  Distortion 


The  simulation  used  in  the  2003  runs  was  adequate  to  the  task  of  determining  the  results 
described  in  the  previous  section:  scientists  at  Thomas  Jefferson  National  Accelerator  Fa¬ 
cility  (JLAB),  where  the  proposed  laser  was  being  developed,  made  use  of  the  information 
provided,  and  the  consideration  of  shorter  Rayleigh  length  configurations  was  taken  up  as  a 
serious  endeavor.4 

For  the  following  year’s  runs,  a  number  of  improvements  were  made  to  the  simulation  code. 
A  faster  fast  Fourier  transform  (FFT)  algorithm5  was  employed  resulting  in  a  10- fold  speed 
improvement.  A  more  accurate  propagation  method  involving  next- nearest  neighbors  was 
utilized,  and  an  expanding  coordinate  system  that  follows  the  rapidly  diffracting  wavefront 
made  for  greatly  improved  efficiency. 

1.  Optical  Mode  Distortion  in  a  Short  Rayleigh  Length  FEL 

A  proposed  high  power  FEL  has  the  parameters  shown  in  Table  3.1.  As  before,  the 
simulation  has  the  ability  to  follow  multiple  and  arbitrary  transverse  optical  modes,  as  they 
interact  with  the  electron  beam  and  bounce  back  and  forth  in  the  optical  cavity,  including 
mirror  transmission  and  edge  losses.6  Output  from  the  improved  simulation  is  shown  in 

Table  3.1:  Parameters  for  a  Proposed  High-Power  FEL 


Parameter 

Symbol 

Value 

electron  beam  energy 

Eb 

80  MeV 

bunch  charge 

q 

400  pC 

bunch  length 

h 

0.15  mm 

peak  current 

i 

800  A 

normalized  emittance 

£n 

3  mm-mrad 

beam  radius 

n 

60  /irn 

undulator  periods 

N 

22 

undulator  wavelength 

^ 0 

2.36  cm 

undulator  length 

L 

52  cm 

peak  magnetic  field 

B 

0.7  T 

undulator  parameter 

K 

1 

optical  cavity  length 

S 

18  m 

Rayleigh  length 

z0 

6  cm 

resonator  quality  factor 

Qn 

4 

optical  wavelength 

A 

1  /irn 

4 At  this  writing,  an  experiment  is  in  place  at  JLAB  involving  a  resonator  configuration  utilizing  a  Rayleigh 
length  on  the  order  of  z  =  0.2. 

5M.  Frigo  and  S.  G.  Johnson,  “FFTW:  An  Adaptive  Software  Architecture  for  the  FFT”,  Proceedings, 
1998  International  Conference  on  Acoustics  and  Signal  Processing,  p.  1381,  1998. 

6  J.  Blau,  et  al,  “Optical  Mode  Distortion  in  a  Short  Rayleigh  Length  Free  Electron  Laser”,  Proceedings, 
26th  International  Free  Electron  Laser  Conference,  Trieste,  Italy,  September,  2004. 
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Figure  3.7.  The  displayed  data  (top  center  of  the  figure)  clearly  exhibits  an  other-than- 
Gaussian  wavefront  for  the  particular  parameters  of  that  simulation  run. 


rn.  HAVEIHOMT  EVOLUTION 


Figure  3.7:  Output  from  a  2004  simulation.  New  visualizations  of  the  optical 
mode  wavefront  accompanied  improved  speed  and  diagnostics. 
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2.  Results 


When  the  wavefront  pictures  provided  by  the  simulation  were  included  in  the  plots  of  the 
parameters  being  varied  to  find  optimum  extraction  (Figures  3.8  and  3.9),  it  became  clear 
that  under  certain  configurations,7  a  FEL  operates  in  higher-order  modes.  It  is  interest- 
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Figure  3.8:  Plot  of  extraction  g  as  a  function  of  dimensionless  Rayleigh  length 
So  and  wavefront  pictures  generated  by  the  simulation. 


Dimensionless  Mirror  Separation,  xmjr 


Figure  3.9:  Plot  of  extraction  g  as  a  function  of  dimensionless  mirror  separa¬ 
tion  rm ir  and  wavefront  pictures  generated  by  the  simulation. 


ing  to  note  that  while  varying  the  Rayleigh  length  Zo  (Figure  3.8)  and  mirror  separation 
rmjr  (Figure  3.9)  evoke  a  relatively  flat  response  from  the  FEL  in  terms  of  extraction,  the 
new  boundary  conditions  imposed  by  these  changes  appear  to  lead  to  a  change  in  modal 
composition  of  the  beam. 

7For  example,  under  certain  Rayleigh  lengths  or  distances  between  resonator  mirrors.  Remember  that 
we  said  the  FEL  automatically  finds  the  optimum  vq,  so  each  data  point  in  Figures  3.8  and  3.9  represents 
an  operating  FEL,  with  associated  pictures  indicating  the  shape  of  its  optical  wavefront. 
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C.  The  Need  for  Higher-Order  Mode  Analysis 


The  radially  symmetric,  generally  papilloid  shape  of  the  higher-order  mode  wavefronts 
observed  in  the  2003  and  2004  simulation  results  is  suggestive  of  the  geometry  of  the  FEL 
resonator  system  and  the  boundary  conditions  imposed  by  an  optical  wave  propagating  along 
and  around  a  coaxial  column  of  current  (see  Figure  3.10,  below).  Under  the  assumption  that 


Figure  3.10:  Wavefront  suggesting  coincident  reinforcement  of  the  optical 
mode  by  the  electron  beam.  This  is  the  general  shape  we  see  in  many  FEL 
simulations,  which  display  the  beam  structure  at  the  high  reflector  end  of  the 
resonator.  The  appearance  of  such  structure  demands  a  detailed  analysis  of 
the  modal  content  of  the  outcoupled  beam. 


the  fundamental  Gaussian  mode  is  the  ideal  mode  both  for  lasing  efficiency  and  atmospheric 
propagation,  a  simple  index  of  beam  quality8  would  be  sufficient.  The  observation  from 
simulations  that  higher-order  modes  occur  within  the  resonator  demands  a  more  robust 
characterization  of  beam  quality.  In  other  words,  we  wish  to  know  more  precisely  which 
modes  comprise  our  beam. 


8One  such  measure  of  beam  quality  frequently  encountered  in  the  literature  is  M2,  given  by 


where  d  is  the  observed  beam  diameter  and 


is  the  diameter  of  a  theoretical  diffraction-limited  Gaussian  beam  arising  from  the  same  initial  conditions. 
For  more  on  M 2  see  Safnett,  and  SlEGMAN  1993. 
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Chapter  IV. 


Higher-Order  Modes 


An  approach  to  the  resolution  of  a  given  laser  beam  into  its  component  modes  is  a  process 
analogous  to  a  Fourier  decomposition  of  a  function.  In  the  same  way  that  an  arbitrary 
function1  may  be  represented  by  a  series  of  sines  and/or  cosines,  an  arbitrary  laser  beam 
may  be  represented  by  a  sum  of  orthogonal  modes.  Identifying  the  appropriate  complete, 
orthogonal  basis  set  is  the  challenge,  and  is  the  topic  of  this  chapter. 


A.  Hermite-Gaussian  Beams 


The  trial  solution  we  used  to  solve  the  parabolic  wave  equation  in  §1.  can  be  further 
generalized  to  allow  for  higher-order  modes.  We  begin  by  noting  that  (2.45)  can  be  expressed 
as  a  product  of  identically  dependent  functions  of  each  transverse  coordinate,  that  is 

S(x,y,T)  =  S0  —  e-^2/^2e^^/(A2+T2)e-iarctan(r/5o) 

W 

—  £Q^e-x2/w2  e~y2/w2  eix2T/{zZ+T2)eiy‘2T/{z%+T2)e-ia,rctan(T/z0) 

W 


arctan(r/zo) j2  ^—x2  j  w2  ^ix2 r / (zq+t2) 


X 


arctan  (t  /  z0) /2  ^-y2 /w2 


eiy2T/(z2+r2) 


=  Sx(x,r)Sy(y,r). 


(4.1) 


As  such,  we  can  simplify  the  process  of  finding  higher-order  solutions  by  examining  only  one 
transverse  dimension  at  a  time.  Writing  the  wave  equation  as 


i  d2  d_\ 
4  dx2  +  8  t  ) 


£x 


0, 


(4.2) 


1  The  function  must,  of  course,  meet  the  Dirichlet  conditions:  on  any  finite  interval,  it  must  possess  one 
minimum,  one  maximum,  and  no  discontinuities. 
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and  modifying  our  trial  solution  to  give  the  complex  amplitude  function  p  an  explicit  trans¬ 
verse  dependence, 

Sx(x,r)  =  £0p{x,T)e~x2/q{T\ 


we  begin  taking  derivatives: 

d2£ 

dx2 


d  2p  4  „  dp 


-Xx:  +  ~zX 


d£ 

~&r 


dx 2  q  dx  \q 

dP  ,  1  d<l  ~2„  \  c  ,-x2/q 


p 


£oe 


-x2/q 


tut  +  —  -j-x  p  £0e 

dx  qz  dr 


to  give 


ijpp+i_idp_  u2  ±  dp  p  dq  ,  iVs  =  0 

4  dx 2  q  dx  q2  2 q  dx  q2  dA 


Dividing  through  by  —  \£oe  x2^q  and  collecting  terms  in  spatial  derivatives  of  p  gives 


d2p  4  _  dp 
dx2  q  dx 


,  ,  M  \ 

1  +  l-r-  ~^X 


dr  J  q 


p  +  4^  =  0. 
dr 


It  is  the  similarity  of  (4.7)  to  the  well  known  Hermite  differential  equation 

d2Hm  dHm 

—  2  u — ; - 1-  2  mHm  =  0 


(4.3) 

(4.4) 

(4.5) 

(4.6) 

(4.7) 

(4.8) 


du 2  du 

that,  after  some  algebra,  leads  to  the  following  expression  for  the  Hermite- Gaussian  beam:2 

V^X  |  —T.^/Zn  Q\ 


£(x,T)=£0^Hm 

w 


e~X2  /  Z0  ei<t>m0(x,T) 


W 


or  in  three  dimensions 


w0 


£(x,y,r)  =  £0—Hm 
w 


'V2: 


X 


w 


Hr, 


'V2  y" 


w 


e-r2/zoei4>mn{r,T) 


(4.10) 


where 

~9 

r t  t 

<t>mn\r,  r )  =  — -  -  (m  +  n  +  1)  arctan  — ,  (4.11) 

Zo  +  T  zo 

and  Hm  and  Hn  are  mth-  and  nth-order  solutions  to  (4.8).  These  so-called  Hermite-Gaussian 
solutions  to  the  dimensionless  parabolic  wave  equation  constitute  a  complete  set  of  orthog¬ 
onal  modes  with  which  one  may  construct  an  arbitrary  laser  beam  wavefront. 

The  properties  of  Hermite  polynomials  are  well  documented.3  The  first  ten  polynomials 
are  listed  in  Table  4.1.  Some  examples  of  the  effect  of  the  functions  Hm  and  Hn  on  the  beam 
wavefront  is  represented  in  Figure  4.1.  These  “burn  patterns”  illustrate  the  intensity  £2  of 
the  beam  in  the  transverse  plane,  but  display  nothing  about  the  phase!  To  fully  characterize 
a  beam,  one  must  know  both  its  amplitude  and  phase  at  every  point,  a  challenge  in  a  real 
experiment,  but  a  matter  of  course  in  a  simulation. 

2The  (dimensional)  functional  forms  of  £(x,y,T)  as  well  as  a  derivation  of  the  higher-order  solution  are 
available  in  A.  E.  Siegman,  Lasers,  University  Science  Books,  Sausalito,  1986.  We  include  the  dimensionless 
forms  of  the  r-dependent  expressions  for  completeness. 

3See,  for  example,  pp.  367-9,  V.  Barcilon,  “Special  Functions”  in  C.  E.  Pearson,  Ed.,  Handbook  of  Applied 
Mathematics,  Van  Nostrand  Reinhold,  New  York,  1983. 
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Table  4.1:  Hermite  Polynomials  -  Solutions  to 


+  2mHm  =  0 


d  Hm  0  d  l  l  h 


Ho{u) 

hm 

H2(u) 

HM 

HM 

HM 

HM 

HM 

HM 

HM 


1 

2  u 

4m2  —  2 
8m3  —  12  u 
16m4  -  48m2  +  12 
32m5  —  160m3  +  120m 
64m6  -  480m4  +  720m2  -  120 
128m7  -  1344m5  +  3360m3  -  1680m 
256m8  -  3584m6  +  13440m4  -  13440m2  +  1680 
512m9  -  9216m7  +  48384m5  -  80640m3  +  30240m 


Figure  4.1:  Some  Hermite-Gaussian  modes.  The  transverse  electromagnetic 
mode  (TEM)  subscripts  correspond  to  the  values  of  m  and  n  in  the  polynomials 
Hm  and  Hn . 


B.  Laguerre-Gaussian  Beams 


In  cylindrical  coordinates,  the  dimensionless  transverse  Laplacian  is 


ld_  (,d\ 

r  dr  \  dr  )  f 2  d62  ’ 


(4.12) 


and  (2.17)  becomes 

IdS  f  dS\  1  d2S  dS 

- r—  +  —  — H - 

r  dr  \  dr  J  r2  dd2  dr 

We  might  guess  that  this  coordinate  system  is  more  appropriate  to  the  resonator  geometry, 
whose  eigenmodes  are  certain  to  reflect  the  circular  shape  of  the  mirrors.  Following  the  same 
procedure  as  in  the  previous  section,  a  generalized  trial  solution  to  (4.12)  yields  a  differential 
equation  similar  to  Laguerre’s  differential  equation 


m- 


d2Lp 
du 2 


+  {m+  1  -  r)  +  pLp  =  0. 
du 


(4.14) 
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Like  the  Hermite  polynomials,  the  polynomials  that  solve  (4.14)  also  constitute  a  complete 
set  of  orthogonal  functions,  the  first  ten  of  which  are  given  in  Table  4.2. 


Table  4.2:  Laguerre  Polynomials  -  Solutions  to  u 


d2Lp 

du2 


+  (1 


\  dLp 


pLp  —  0 


L0(u ) 
Li(u) 
L2{u ) 
L3(u) 

La[u) 

L5(u) 

Lq{u) 

L7(u) 

L8(u) 

Lg(u) 


1 

1  —  u 


1  —  2  u  +  \u2 
1  —  3  u  +  f  «2  —  |w3 
1  —  4-u  +  3  u2  —  |  u3  +  ^u4 

1  -  5u  +  5w2  -  | u3  +  £ u 4  -  j^u5 
1  _  6u  +  fu2  -  +  |m4  _  i_us  +  _i_Me 

l-7u+fu2-fu3  +  lu4-^5  +  ^6-^7 
1  -  8U  +  Uu2  ~  f  “3  +  if  “4  -  J5U 5  +  T80m6  -  <ku7  +  lok^8 
1  -  9u  +  18“2  -  14«3  +  f uA  -  1“5  +  ~  Tio^7  +  4 -  364om9 


The  Laguerre- Gaussian  beam,  given  by 


T(r,0,r) 


C  ^0 

C-0~ 

w 


e-r2/u>2 eic/>™(r,e,T)  ^ 


(4.15) 


where 

~2 

0™(r,  0,  r)  =  „r  T  -  —  (2p  +  m  +  1)  arctan  —  +  m6,  (4.16) 

p  -5  + r  zo 

is  a  dimensionless  solution  to  the  parabolic  wave  equation  in  cylindrical  coordinates.4  Beam 
propagation  is  shown  for  the  first  ten  Laguerre-Gaussian  modes  in  Figure  4.2.  Evidently,  the 
wavefronts  retain  their  shape.  Examining  (4.15),  we  see  there  is  no  change  in  the  number  of 
zeros  as  r  increases:  there  is  no  “crossing”  of  local  maxima. 

4The  subscripts  ?n  and  n  are  standard  in  the  literature  for  the  Hermite  polynomials.  Similarly,  m  and  p 
are  often  used  for  the  Laguerre  indices.  At  this  point,  we  essentially  leave  behind  our  discussion  of  Hermite- 
Gaussian  beams,  so  from  now  on  m  will  denote  the  angular  mode  of  the  Laguerre-Gaussian  solution. 
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Figure  4.2:  The  first  ten  Laguerre-Gaussian  modes.  Notice  that  L §  is  simply 
a  Gaussian  beam.  Notice  also  that  all  modes  retain  their  general  shape — in  par¬ 
ticular  the  number  of  zeros — as  they  propagate  through  a  sourceless  medium. 


The  orthogonality  of  the  generalized  Laguerre  polynomials  is  weighted  by  ume  11 
normalization5 

[°°  ume-uL;  (u)L^(u)du  =  r(g  +  7  +  1)<W 
Jo  Q- 


with  the 
(4.17) 


where 


r(a)  = 


e  uva  1du. 


(4.18) 


Jo 

We  shall  only  consider  those  Laguerre  modes  with  no  azimuthal  dependence,  that  is  m  =  0. 
For  integers  (which  q  +  1  surely  is) 


r(n)  =  (n  —  1)!, 


(4.19) 


and  so,  for  our  purposes,  we  shall  write  the  orthogonality  condition  as 

e~uLp(u)Lq(u)du  =  5pq.  (4.20) 

Now,  hireling  the  Laguerre-Gaussian  (LG)  “spectrum”  of  an  arbitrary  initial  wavefront 
E(r,  0)  with  azimuthal  symmetry  requires  us  to  calculate  the  coefficients  cp  in  the  series 

E(f,  0)  =  cp£p(?,  0)  (4.21) 

p 

5See  pp.  533-4  in  D.  Zwillinger,  CRC  Standard  Mathematical  Tables  and  Formulae ,  Chapman  &  Hall, 
New  York,  2003. 


43 


where 


£P(r,  0)  Sq Lp  _  I  e 


2  f5 


-r+io 


(4.22) 


solves  (4.13)  at  r  =  0.  Notice  that  the  factor  e in  (4.15)  becomes  unity  since 

y X  0  Q 

0°(r,  0,  0)  =  7\ - (2 p  +  0  +  1)  arctan - h  0  x  0 

py  ’  zg  +  0  v  y  ’  z0 

=  0. 

Expressing  onr  arbitrary  wavefront  as  a  function  of 


Q  = 


2  r2 


and  noting  that6 


w- 


_  2 r2  _  2 r2 

Qo  ~  o  ~  5 

Wq  Zq 


we  multiply  both  sides  of  (4.21)  by  an  LG  mode  £q(go)  and  integrate  over  g0: 

POO  POO 

/  E(g0)£q(go)dg0  =  /  ^jcp£p(g0)£q(go)dg0 

Jo  Jo  p 

POO 

=  ^p  /  e~Q°Lp(Qo)Lq(go)dg0 

p  J° 

=  £'o  CpSpq 


(4.23) 


(4.24) 


(4.25) 


-  £2c 


6For  example,  using  this  new  coordinate  the  Gaussian  beam  becomes 

£0e~f2/~z°  ->  £0e~eo/2. 


(4.26) 
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where  we  have  used  (4.20).  An  arbitrary  coefficient  is  then' 


Cq  £2  /  E(,Qo)£q(,Qo)d'Qo 


'0  Jo 

4 


£(Ao 


fe~r2/Z0Lq  (  )  E(r)dr. 

Zq 


(4.27) 


Equation  (4.27)  is  limited  to  r  =  0.  Even  as  we  assert  the  usefulness  of  a  modal  analysis 
of  the  initial  wavefront,  we  wish  to  analyze  our  beam  at  any  point  along  its  propagation, 
that  is  for  any  r.  Beginning  once  again  with  (4.21),  this  time  we  make  no  assumptions  about 
the  value  of  r  and  thus  w(r)  and  the  r-dependent  phase  factor  el<f>p’^r,e,T'>: 


E{q,t)  =  ^cpEp(q,t) 
p 


(4.28) 


Implicit  in  this  generalization  is  the  possibility  of  a  complex  cp(r).  Here  is  why.  Suppressing 
the  arguments  g  and  r,  let  us  multiply  both  sides  of  (4.28)  by  arbitrary  mode  E*  and 


'We  can  normalize  the  quantity  JAcp  to  unity  by  setting  E(f)  to  the  fundamental  and  calculating 


where  we  have  used 


Co  = 


Eo  So  Jo 

4  r 


£0z0 

4 


re-f2/Soe-f2/iodr 


re~2f2/'Z0dr 


'o 


_ Zo 

SqZo  4 

1 

fo’ 


T[(n+  1  )/p\ 
pa(n+1)/p 


with  n  =  1,  a  =  2/zq,  and  p  =  2  (see  p.  451,  Zwillinger).  Having  thus  identified 

Co 

and  returning  to  our  standard  LG  index  p,  a  normalized  Laguerre-Gaussian  spectral  coefficient  is  then 

'  2r2' 


/2o  j 


Zo 


E(f)dr, 


with  the  carat  indicating  normalization.  While  this  is  a  useful  result,  we  ultimately  intend  to  approach  the 
calculation  of  LG-spectrum  coefficients  numerically,  and  a  self-consistent  normalization  will  be  included  in 
the  computer  program.  The  above  is  included  for  completeness. 
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integrate: 


poo  POO 

/  E£*dg  =  /  ^2 CpSpSqdg 

Jo  Jo 


£oY Zcp  /  e*(0P  ^)e  °Lp(g)Lq(g)dg 

P  Jo 

POO 

^^-,1^°)  /  e-eLp(Q)Lq(Q)dg 

P 

c2  2i(p—q)  arctan(r/i0)  X 

°0  /  y  CPC  upq 


-  £2c 
c'0Lqi 

and  so,  switching  indices  from  q  to  our  standard  p, 


cp  =  p^  E£*dg 


-9  c  2 


or  explicitly 


/*uu 

c  =  ^Lei(2p+l)arctan(r/«)  /  e-3e/2 L  (g)E(g,T)dg 

£qw  Jo 


_±ei{2P+l),rctMr/z0)  /  fe~^  Ivd  L  '  E(f,r)df, 

£0w  J  o  \r/ 


where  once  again 


3  =  1  +  ^- 


(4.29) 


(4.30) 


(4.31) 


(4.32) 


We  will  attempt  to  solve  (4.31)  numerically,  for  two  reasons:  (1)  while  we  may  know  the 
amplitude  and  phase  of  a  simulated  FEL  beam  at  every  (finite)  transverse  point,  we  will  not 
generally  have  an  analytic  expression  for  the  entire  wavefront;  and  (2)  even  if  we  did  know 
the  functional  form  of  E(r,r),  the  integral  in  (4.31)  is  difficult  to  perform  for  all  but  the 
most  ideal  beams. 

For  a  numerical  analysis,  we  must  first  get  the  integral  into  a  form  the  computer  can 
understand.  Letting  the  subscripts  r  and  i  indicate  the  real  and  imaginary  parts,  we  rewrite 
(4.30)  as 

1 

cP  =  j2  £vEdQ 


c2  /  'i£p,i)  (Er  iEi)  dg 


£20  Jo 


c2  /  (£p,rEr  i£p,rEi  i£p,iEr  -l-  £p^EjJ  dg 


£20  Jo 


c2  /  (£p,rEr  +  £p,iEi)  dg  +  2  /  ( £p,rEi  £p^Er^  dg 


£20  Jo 


£20  Jo 


(4.33) 


We  may  now  split  cp  into  its  real  and  imaginary  parts,  each  of  which  may  be  integrated 
separately: 


£p,r  — 

i 

/»oo 

/  (£p,rEr  +  SppEi) 

J  0 

i  dg, 

(4.34a) 

hi 

i 

°0  ■ 

/»oo 

/  (£  E  —  £  FT 

Jo 

1  dg. 

(4.34b) 

The  contribution  to  the  LG  spectrum  of  an  arbitrary  beam  by  an  individual  LG  mode  shall 
be  given  by  the  modulus  of  that  mode’s  coefficient, 


CP  \Cp\  \J  Cp,r  “b  Cp,H 

and  will  be  normalized  to  the  sum  of  all  basis  set  coefficients: 


(4.36) 


C.  Numerical  Propagation  of  Higher-Order  Modes 

The  diagnostic  tool  to  be  included  in  free  electron  laser  simulations  numerically  integrates 
(4.34)  and  reports  cv  for  each  of  the  first  ten  LG  modes  via  (4.35).  Figure  4.3  is  a  sample 
output  from  bm.c,  a  self-contained  code  that  resolves  an  input  laser  beam  wavefront  into 
component  LG  modes  and  reports  the  normalized  LG  coefficients  as  a  bar  chart  beneath 
the  wavefront  picture  both  before  and  after  propagation.  In  this  case,  the  integration  has 
correctly  identified  the  initial  wavefront  as  an  exclusively  LG-2  beam:  in  the  bar  chart  below 
the  initial  wavefront  intensity  plot,  only  the  bar  in  the  third  position  is  present  and  goes 
to  the  top  of  the  chart.  To  within  a  few  percent  error,  the  program  has  also  determined 
the  final  wavefront  to  be  LG-2,  as  we  expect.  The  text  in  the  middle  of  the  output  reports 
parameters  used  for  that  simulation  as  well  as  a  “fractional  energy  on  target”  for  a  circular 
target  of  fixed  radius.  This  is  discussed  in  the  next  section.  The  basis  sets  used  to  analyze 
the  input  beam  both  before  and  after  propagation  are  evaluated  analytically  and  displayed 
at  the  bottom  of  the  output  for  reference. 

The  diagnostic  is  able  to  pick  out  multiple  modes,  as  in  Figure  4.4.  In  this  case,  we  see  for 
the  first  time  the  “crossing”  of  zeros  in  the  propagation  of  the  laser  beam.  Remember  that 
“pure”  LG  modes  do  not  exhibit  this  behavior:  a  given  Laguerre-Gaussian  mode  possesses 
the  same  number  of  zeros  throughout  its  propagation  (see  again  Figure  4.2).  Single  modes 
also  always  exhibit  a  global  maximum  at  f  —  0,  not  always  observed  in  mixed  modes. 
This  distinction  between  single-  and  mixed-mode  higher-order  beams  immediately  suggests 
a  useful  application  for  our  tool,  and  it  leads  us  to  consider  one  of  the  difficulties  in  the 
experimental  analysis  of  laser  beam  modes. 

It  is  not  possible  to  determine  the  modal  composition  of  a  beam  (function)  if  only  its 
intensity  (square)  is  known.  As  a  simple  illustration,  let  us  consider  the  curves  in  Figure  4.5. 
It  is  not  possible  to  determine  whether  the  curve  on  the  right  was  arrived  at  by  squaring  the 


(4.35) 
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Time  Steps  (ts)  =  20000  (13888  required) 
Wndow  Size  (W)  =  18.000000 
"Vfeist  Location  (tauw)  =  0.000000 
Rayleigh  Length  (zO)  =  0.500000 
Number  of  Bins  (n)  =300 
Target  Radius  (wm)  =  0.450000 
Fractional  Energy  on  Target  =  0.107658 


C=0.003472 
error =0.0 3 03 66 


Figure  4.3:  Output  from  the  diagnostic  tool  bm.c.  The  basis  set  before  and 
after  propagation  (evaluated  analytically)  is  included  at  the  bottom  of  the 
display  for  reference. 


top  left  or  bottom  left  function.  An  attempt  to  determine  the  modal  composition  of  a  laser 
beam  based  on  intensity  information  alone  encounters  the  same  problem,  ft  is  precisely  the 
phase  information  that  is  lost  in  the  process  of  squaring. 

Consequently,  burn  patterns  can  be  misleading.  Figure  4.6  shows  the  propagation  of  the 
same  laser  beam  as  Figure  4.4,  but  takes  the  cross  sectional  amplitude  (or  burn  pattern)  of 
the  final  wavefront  earlier  than  before.  Visually  inspecting  the  initial  and  final  wavefront 
amplitude  plots  (enlarged  in  Figure  4.7),  one  might  be  hard-pressed  to  conclude  (or  even 
guess!)  that  they  were  generated  by  the  same  laser  beam.  Yet  the  values  of  the  LG  coeffi¬ 
cients,  as  reported  in  the  bar  charts  beneath  the  initial  and  final  wavefronts  in  Figure  4.6, 
assure  us  that  this  is  the  case.  Multiple  modes  interfere  because  each  mode  prescribes  a 
unique  time  dependence  via  the  factor 


eiMr,r)  =  exp 


r2r 


r 


Zn  +  A 


L  ^0 


—  i  (2p  +  1)  arctan  — 

z0 


(4.37) 


in  (4.15).  This  evolving  phase- front  results  in  an  evolving  amplitude-front  as  modes  cyclically 
cancel  and  reinforce  each  other,  as  evident  in  the  propagation  plot. 

The  diagnostic  tool  performs  well  for  even  more  complicated  mixtures  of  higher-order 
Laguerre-Gaussian  modes.  In  Figure  4.8,  bm.c  has  correctly  identified  the  3:5:2  mixture  of 
the  LG-0,  -3,  and  -6  modes  comprising  the  simulated  beam.  Figure  4.9  shows  the  analysis 
of  a  beam  comprised  of  many  LG  modes.  Once  again,  in  both  figures  we  observe  the  zero 
crossings  typical  of  mixed-mode  beams,  and  the  difference  in  the  burn  patterns  before  and 
after  propagation. 
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Figure  4.4:  A  beam  composed  of  two  distinct  Laguerre-Gaussian  modes.  The 
slight  error  in  the  bar  chart  pertaining  to  the  LG  coefficients  for  the  final 
wavefront  is  a  result  of  numerical  integration  error. 


Figure  4.5:  When  squared,  the  two  functions  sin  a:  (top  left)  and  sin(a;  +  tt) 
(bottom  left)  lead  to  an  identical  squared  function  (right). 


D.  Power  on  Target 

Knowing  which  modes  comprise  our  laser  beam  allows  us  to  examine  the  relationship 
between  physical  properties  of  the  beam  and  its  composition.  An  example  of  this  is  the 
important  consideration  of  power  on  target:  how  much  of  the  beam’s  initial  energy  is  focused 
onto  a  target  at  the  propagation  distance. 

Choosing  a  circular  target  of  an  arbitrary  radius  such  that  almost  all  (93.5%)  of  the 
outcoupled  energy  is  delivered  to  its  surface  when  the  beam  is  in  the  fundamental,  Figure  4.10 
plots  the  energy  delivered  to  the  same  target  by  beams  in  unmixed,  higher-order  Laguerre- 
Gaussian  modes. 

Evidently,  by  themselves  higher-order  modes  are  not  good  choices  for  delivering  power  to 
a  target  at  the  center  of  the  beam.  As  p  increases,  more  and  more  energy  is  relegated  to  the 
“shoulders”  of  the  central  field  (see  Figure  4.11).  Once  in  this  outer  region,  light  cannot  cross 
back  into  the  center  without  a  change  in  the  number  of  local  field  maxima — not  observed  in 
pure  LG  beams. 

Certainly,  real  lasers  are  unlikely  to  be  comprised  solely  of  single  higher-order  Laguerre- 
Gaussian  modes,  or  even  predictable  mixtures.  Figure  4.12  shows  three  theoretical  lasers: 
one  a  Gaussian,  one  an  even  mixture  of  the  Gaussian  and  LG-1  modes,  and  one  a  beam  with 
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Figure  4.6:  Each  LG  mode  carries  a  different  phase  factor,  causing  periodic 
interference  when  multiple  modes  are  combined. 


Figure  4.7:  A  combination  of  LG  modes  before  (left)  and  after  (right)  prop¬ 
agation  over  a  finite  distance  (enlarged  from  Figure  4.6).  While  the  two  wave- 
fronts  are  not  necessarily  identifiable  by  eye  as  manifestations  of  the  same  laser 
beam,  the  diagnostic  assures  us  this  is  so. 


a  “flat-top”  cross  section.  The  percentage  of  energy  received  by  the  centrally  located  target 
is  shown. 

What  the  propagation  of  mixed  modes  such  as  the  flat-top  beam  teaches  us  is  that  the 
inclusion  of  higher-order  in  the  laser  does  not  necessarily  spell  doom  for  the  propagation  of 
power  to  the  target  area.  Certainly  the  trend  in  Figure  4.10  is  not  towards  larger  central  held 
amplitude  for  higher  values  of  p.  But  mixed  modes,  as  Figures  4.6,  4.8,  and  4.9  demonstrate, 
have  the  potential  to  get  some  of  that  power  from  the  shoulders  back  into  the  center  of  the 
beam.  In  fact,  we  may  conclude  that  any  non-fundamental  beam  delivering  more  power 
to  our  test  target  than  the  L\-Gaussian  is  by  definition  comprised  of  multiple  higher-order 
modes.  Particularly  for  use  as  a  weapon,  if  a  laser  beam  is  not  in  the  fundamental,  it  may 
benefit  from  the  mixing  of  higher-order  modes. 
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Figure  4.8:  A  beam  comprised  of  three  higher-order  modes.  The  recipe  is 
three  parts  LG-0  (the  fundamental  Gaussian),  five  parts  LG-3,  and  two  parts 
LG-6.  Once  again,  the  burn  patterns  before  and  after  propagation  are  difficult 
to  reconcile  visually. 
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r-n  i  i  l~kr 


Figure  4.9:  A  beam  comprised  of  many  higher-order  modes. 


p 


Figure  4.10:  EO  and  El  are  the  energy  contained  in  the  beam  and  the  en¬ 
ergy  delivered  to  the  target  area,  respectively.  For  unmixed  Laguerre-Gaussian 
modes,  laser  power  on  target  drops  off  significantly  for  p  >  0. 
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Energy  in  the  “shoulder”  falls  outside  the  target  area. 


Figure  4.11:  A  target  just  large  enough  to  receive  nearly  100%  of  the  energy 
in  the  fundamental  will  generally  not  receive  energy  from  pure  higher-order 
modes  beyond  the  first  zero-point. 
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Figure  4.12:  The  propagation  of  (a)  a  Gaussian  beam;  (b)  a  beam  comprised 
of  equal  parts  L0-  and  L\ -Gaussian  modes;  (c)  an  initially  “flat-top”  beam, 
defined  by  a  step  function  at  f  =  wo  as  determined  by  the  fundamental  Gaus¬ 
sian  in  (a).  Notice  once  again  the  crossing  of  local  maxima  in  (c),  not  observed 
for  unmixed  Laguerre-Gaussian  modes.  It  is  this  crossing  that  has  allowed  the 
beam  to  deliver  more  energy  to  the  center  of  the  target. 
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E.  Conclusion 


Resolving  arbitrary  optical  modes  into  higher-order  mode  spectra  provides  a  useful  tool 
for  the  analysis  of  FEL  simulations  and  a  self-consistent  index  of  beam  quality.  Future  work 
might  include  retaining  the  azimuthal  index  m  in 


£(f,e,T) 


C  U>0 

to  — 
w 


e-r2/w2ei<t>™(r,e,r ) 


where 


~9 

r  r  t 

-5 - -  —  (2 p  +  m  +  1)  arctan - 1-  m9, 

Zq  +  T2  Z0 


during  the  analysis — such  a  basis  set  would  be  capable  of  generating  any  beam,  even  those 
without  radial  symmetry.  Free  electron  laser  simulations  have  also  demonstrated  the  poten¬ 
tial  usefulness  of  the  Hermite-Gaussian  modes  as  an  appropriate  basis  set  when  studying 
off-axis  electron  beams  and/or  resonator  cavity  mirror-tilt. 
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Appendix:  The  Parabolic  Wave 
Equation 


We  present  the  parabolic  wave  equation  in  both  dimensionless  and  dimensional  form. 
Cartesian  coordinates. 

The  dimensionless  parabolic  wave  equation  is 


where 


88  i  2  ~ 

Tr  =  4V^’ 


z 

T  =  V 


d2  d2 

Vi  = 


dx2  dy2  ’ 

for  a  propagation  distance  L.  Using  the  dimensionless  coordinates 


x  = 


LX 


x , 


we  can  transform  to  the  dimensional  paraxial  wave  equation  as  follows: 


and  finally  the  familiar 
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Cylindrical  coordinates. 


55 


In  cylindrical  coordinates  we  use 


V2  ^  + 1£- 
_L  r  dr  \  dr  J  r2  d62  ’ 

where  r  =  r /wl  and  6  =  6 ,  and  we  have 

i  Id  fdS\  1  da 2 
4  f  dr  \  dr  J  r2  d62 
iZX  ri  d  f  dS\  1  dS2' 

47t  r  dr  \  dr  J  +  r2  d62 

i  ri  d  (  dS\  1  d£2 1 
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and  finally 
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